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Abstract 



We consider a branching particle system where each particle moves as an 

" ' independent Brownian motion and breeds at a rate proportional to its distance 

from the origin raised to the power p, for p G [0, 2). The asymptotic behaviour of 

I— I the right-most particle for this system is already known; in this article we give 

P^ large deviations probabilities for particles following "difficult" paths, growth 

Ph rates along "easy" paths, the total population growth rate, and we derive the 

cl optimal paths which particles must follow to achieve this growth rate. 
-(— > 

a 

S 1 Introduction and heuristics 

(N 1.1 The model 
> 

CO We study a branching Brownian motion (BBM) in an inhomogeneous breeding po- 

^~r tential on R. Fix /? > 0, p G [0,2), and a random variable A, which takes values 

^^ in {2,3, . . .}, satisfying E[yllog^] < oo. We initialise our branching process with a 

j^ single particle at the origin. Each particle u, once born, moves as a Brownian mo- 

<^ tion, independently of all other particles in the population. Each particle u alive at 

^^ time T dies with instantaneous rate (3\Xu{T)\p , where X„(T) is the spatial position 

• • of particle u (or of its ancestor) at time T. Upon death, a particle u is replaced 

.5^ by a random number A^ of offspring in the same spatial position, where each A^ 

^ is an independent copy of A. We define m := K[A] — 1, the average increase in the 

^ population size at each branching event. We denote by N{T) the set of particles 
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alive at time T. We let P represent the probability law, and E the corresponding 
expectation, of this BBM. 

The case p = 2 is critical for this BBM: if the breeding rate were instead /3| • I'' 
for p > 2, it is known from Ito and McKean |21j that the population explodes in 
finite time, almost surely. For p = 2, the expected number of particles explodes in 
finite time, but the population remains finite, almost surely, for all time. 

Branching Brownian motions are closely associated with certain partial differen- 
tial equations. In particular, for the above BBM model, the McKean representation 
tells us that 



solves the equation 



v{T,x):=e[ H f{Xu{T)) 
dv 1 d V 

a^ = 2ij^ + W(GW-.) (1) 

with the initial condition v{0,x) = f{x), where G{s) := E(s ) is the generating 
function of the offspring distribution A. In the case of constant branching rate {p = 
0), this is known as the Fisher-Kolmogorov-Piscounov-Petrovski (FKPP) reaction- 
diffusion equation. 

An object of fundamental importance in the study of branching diffusions is 
the right-most particle, defined as Rt '■= max„g7v(2i) X^(r). Standard BBM, with 
binary branching at a constant rate (that is, p = and G{s) = s"^), has been much 
studied. In this case, it is well known that the linear asymptotic limx'_>oo Rt/T = 
\/2f3 holds almost surely. The distribution function of the right most particle position 
solves the FKPP equation with Heaviside initial conditions, and it is known that 
¥{Rt > m(T) + x) — )■ w{x) where w is a travelling-wave solution of M and m{T) 
is the median for the rightmost particle position at time T. Sub-linear terms for 
the asymptotic behaviour of m{T) = y/2j3T — 3/(2y^2^) logT + 0(1) were found by 
Bramson ^ and [7j . See also the recent shorter probabilistic proofs by Roberts ^26j , 
and corresponding results for branching random walk by Aidekon [Ij and Hu and 
Shi [2D]. For approaches using partial differential equation theory, see the recent 
short proof by Hamel et al. [T2] and an impressive higher order expansion due to 
Van Saarloos [28j . Detailed studies of the paths followed by the right-most particles 
have been carried out by Arguin et al. [31 Hj , and by Aidekon et al. [2j . 

For p G (0, 2), right most particle speeds much faster than linear occur and Harris 
and Harris [T6] found an asymptotic for Rt using probabilistic techniques involving 
additive martingales and changes of measure. 

Theorem 1 (Harris, Harris [16]): 
Forpe [0,2), 



T^oo 2^2^ V 2 



1 



l.m^=m2-p)=V- 



almost surely. 



(Note that, as above, the theorems given in the present paper are all written for 
p G [0, 2). However most of our results were already known in the case p = 0.) 



In this paper we study in more detail the paths followed by particles in the 
BBM. Theorem [T] suggests a rescaling of time and space, and we consider whether 
particles follow paths which, after rescaling, lie in a particular subset of C[0, 1]. 
In section [2| we give large deviations probabilities for particles following "difficult" 
rescaled paths as well as results on the almost sure growth rates for the number of 
particles following any "easy" path. From these results we can derive the growth 
rate of the total number of particles in the BBM, and find the paths which particles 
must follow to realise this growth rate; this involves solving certain path optimisation 
problems subject to integral constraints, the solutions to which are not obvious, but 
nevertheless can be found explicitly and have intuitive probabilistic interpretations. 
A surprising and very significant feature arising from the (unbounded) spatially 
dependent branching rate of this model is the fact that the expected number of 
particles and typical number of particles following paths do not match, even on the 
exponential scale. 

Although this work is the natural sequel to |16j . spatially dependent branching 
rates have not often been studied in detail. See Git et al. [TT] , and Lalley and Sellke 
|221 [23] for a case with bounded breeding potential. Other studies of branching 
processes with time inhomogeneous environments include recent works by Fang and 
Zeitouni [9l[10j, where analogous path optimisation problems also appear. Recent 
developments in the study of spatially inhomogenous versions of the FKPP equa- 
tion from the PDE's perspective include, for example, the periodic environments in 
Hamel et al. |13| . A key technique which is used in |121I13) is to relate the non-linear 
PDE problem to a free boundary linearised PDE one. In fact, such free boundary 
problems are intimately related to the probabilistic constrained path optimisation 



problems (discussed in Section 2.4). We note that, with the exception of [Ml [TT] . 
the above articles are all concerned with bounded environments. 

Unbounded branching rates lead to unusual features and pose considerable tech- 
nical difficulties, much as their corresponding unbounded non-linear differential op- 
erators would. One manifestation of the unbounded branching rates is the position 
of the right-most particle growing faster than linearly in time (as in Theorem [I| ; 
another is the disagreement of expected and typical particle behaviours. 

We start by giving a very rough heuristic explanation for some of our results. The 
technical details, relevant definitions and precise statements of our main results will 
be given in section [2J We have strived to make the heuristics as clear as possible in 
the hope that the reader can gain a good understanding of our main results without 
necessarily having to read the technical details in the rigorous proofs of sections 
[3]|6j The solution of the constrained path optimisation problem and it properties 
are found in sections [7] and [H 

1.2 Heuristics 

Whilst we are mainly interested in the almost sure behaviour of the inhomoge- 
neous branching Brownian motion, this is typically much harder to obtain than the 
expected behaviour. However, we will be able to get a very good intuitive under- 
standing of the almost sure behaviour by carefully considering expectations of the 



number of particles that travel close to given trajectories. 

In a sense that Schilder's Theorem |29] from Large Deviation theory can make 
precise, the small probability that a Brownian motion B manages to stay 'close' to 
some given trajectory F is very roughly given by 

F{B stays 'close' to F during [0,r]) w exp ( -- / F'(s)^ds j . 

where F : [0, T] — )• M is suitably 'nice' with F{0) = and T is some large time. 

Since any particles that are close to trajectory F at time s will give birth to an 
average of m new offspring at a rate close to /3\F(s)\p , the expected total number 
of particles that have stayed close to trajectory F up to time T will very roughly be 
given by 



E[#{u G N{T) : X„ 'close' to F}] ^ exp ( f 



m/3|F(5)r- V(s)2 



ds). (2) 



=)^TA/(i) 






(3) 


is some fixed function, then (2) leads to 




2+p /■! 
f} ~ T2-P / 
Jo 


mPlfisr-^fisf 


ds, 


(4) 



Heuristically, we can think of the number of particles travelling along a 'nice' tra- 
jectory F as behaving roughly like a time-dependent birth-death process (see [19]) 
with a birth rate mP\F[s)\P and a death rate ^F'{s)'^. It is now natural to look for a 
scaling of paths where the birth and death rates are of the same order of magnitude. 
That is, if we consider trajectories of the form 

F(s) 

where T is large and / : [0, 1] — ; 

logE[#{n G N{T) : u follows /}] 

where by 'u follows /' we mean that the rescaled particle path, T^^'(^^P'X„(sT), 
remains within some very small distance of the given rescaled path /(s) for all 
s £ [0, 1]. Essentially, this is Theorem ^ which reveals how the expected number of 
particles varies along a given scaled up path: heuristically, for t G [0, 1] and T large, 

logE[#{u £ N{tT) : u follows / up to time t}] ~ T^K{f,t) (5) 

where 

K{f,t):= f\mf3\fisr-l^f{sf]ds (6) 

for 'nice' functions / : [0, 1] — )• M with /(O) = 0. The functional K{f, t) is of 
fundamental importance for the inhomogeneous BBM and will play a crucial role 
throughout this paper. 

We will show that the expected number of particles which end up near the 
rescaled position z (corresponding to actual position T'^i\^~pi z) grows like the ex- 
pected number of particles following some optimal rescaled path hz ending at z, 

r n 2+p 

logE[#{ii E N{T) : rescaled path of u ends near z}] ~ T^-pK{hz, 1), (7) 



where K{hz,l) = supf {K{f,l) : /(I) = z}. The optimal function hz satisfies 
the corresponding Euler-Lagrange equation h'!, + mph^ = with hz{0) = and 
hz{l) = z. Optimising over z then suggests that the expected total population size 
satisfies 

logE iV(T) ^ T^^ sup K{hz, I) =T^^snpK{f,l). (8) 

/ 

Theorem |4] confirms these heuristics and explicitly identifies the expected total pop- 
ulation growth rate. 

However, a more surprising fact is that these results on the expected number 
of particles are not representative of a typical realisation of the system. Indeed, a 
dominant contribution to the expected number of particles at large time T can come 
from vanishingly rare events where particles go very far away from the origin to take 
advantage of high reproduction rates. Note that, for any t G [0, 1], 

P(some particle follows / up to t) = inf P(some particle follows / up to time s) 

se[o,t] 



< inf E\#{u G N(sT) : u follows / up to s}] , 



se[o,t] 
^exp 



T^ mi K{f,s)\. (9) 



If infggp^i] -ftr(/, s) < 0, then the path / is 'difficult' and it is very unlikely to be 
observed in a typical realisation, even if K{f, 1) > so that the expected number 
of particles alive at time T having followed that path is very large. We are in fact 
able to show that, for a difficult path, the last inequality in ([9]) will be attained up 
to leading order in the exponent, and hence that 'difficult' paths will satisfy 

2+p 

log P (some particle follows / up to t) ~ T'^-v inf K{f, s) < 0. (10) 

se[o,t] 

This probability of presence result is stated rigorously in Theorem |8j Looking for 
the trajectory that travels the furthest without being 'difficult' leads us to guess 
that the right-most particle boundary satisfies mjir^sY — 2^'{s)'^ = with r(0) = 
(in agreement with Theorem [l]). 

On the other hand, if we have not yet had any 'difficult' points along the path, 
we might guess that the almost sure and the expected growth rates will still agree. 
Indeed, Theorem [5] will confirm that, roughly speaking, we almost surely have 

log#{n G N{tT) : u follows / up to t} 



2+p 



^T^^K{f,t) if K(/,s) >Ofor ahsG [0,t]. (11) 

for any t G [0, 1]. Hence, if K(f,s) is always non-negative there will almost surely 
be some particles following / in the large T limit. On the other hand, if K{f,s) 
becomes strictly negative for the first time at some time ^o then it becomes expo- 
nentially unlikely that any particle makes it past this bottleneck; Oq corresponds to 
the extinction time along this 'difficult' path /. 



Finally, we anticipate that the almost sure number of particles which end up 
near the rescaled position z, with \z\ < r(l), will grow like the expected number of 
particles following some optimal path Qz that does not undergo any extinction: 

log#{n G N{T) : rescaled path of u ends near z} ~ T^~p K{gz, 1), (12) 

where K{gz,l) = supf{K{f,l) : /(I) = z,K{f,s) > for ah s G [0, 1]}. The op- 
timal path gz that gives rise to the vast majority of particles turns out to have 
two distinct phases. Initially, it follows the trajectory of the right-most parti- 
cle, thereby gaining optimal potential for future growth without becoming extinct. 
Then, after some optimal intermediate time, it "cashes in" during the second phase, 
switching over to the path that maximises growth, which is the path that satisfies 
g'z + mpg^z~ = with 52(1) = z. We will see that, in addition, the optimal path 
necessarily has a continuous derivative and this property determines the point at 
which one switches from one phase to the other. 

The almost sure total number of particles in the system can now be recovered 
by optimising over z, with 

2+p 2+p 

log|iV(r)| ~r^sup/^(g2,l) = r^sup{K(/,l) ■.K{f,s) > for all s G [0,1]}. 

/ 

In particular, this will reveal that the almost sure population growth rate is strictly 
smaller than the expected population growth rate; the constraint that the paths 
cannot have passed through any extinction times has a significant effect. 

2 Main results 

Fix a set -D C C[0, 1] and t £ [0, 1]. We are interested in the sets 

NT{D,t) := j-u G N{tT) :3f £D with Xu{sT) = T^f{s) Vs G [0,t]} (13) 

for large T. We will typically consider sets of the form D = B{f,e) for a given / 
(the ball of C[0, 1] with centre / and radius eYj in this case NT{D,t) is the set of 
particles alive at time tT whose rescaled paths up to that point have stayed within 
distance e of /. Thus, for a given path / which we keep rescaling in space and time 
according to large T, Nt{D, t) tells us how the population following that path grows 
and shrinks as t varies between and 1. 
Define the class Hi of functions by 

ifi := |/ G C[0, l]:3ge L^[0, 1] with f{t) = j g{s) ds Vt G [0, 1] | , 

and to save on notation set f'[s) := cxd if / G C[0, 1] is not differentiable at the point 
s. Observe that f £ Hi implies /(O) = 0. 



^In this paper, the space C[0, 1] of continuous functions on [0, 1] is always endowed with the La 
topology. 



We can now define K precisely: for / G C[0, 1] and t G [0, 1], 



KU,t)--- 




t r- 1 

otlierwise. 



mfi\f{s)\^--J{s) 



ds if/GFi, 



We use throughout the paper the convention that inf = +oo and sup0 = — oo. 

2.1 Expected population growth 

Our first result is rather straightforward and gives the behaviour of the expectation 
of the number of particles following paths in some set. 

Theorem 2: 

For any closed set D C C[0, 1] and t G [0, 1], 

hmsnp^logE\NTiD,t)\ < supi^(/,t), 

and for any open set A C C[0, 1] and t G [0, 1], 

liminf ^:^logE|iVT(yl,t)| > sup K{f,t) 

Moreover, if we define 

i^E(^):=sup{i^(/,l):/GC[0,l], /(I) = ^}, (14) 

we have the following easy corollary: 

Corohary 3: 

For each e > and z G M, let D^^, := {/ G C[0, 1] : |/(1) - z\ < e}. Then 

lim lim ^logE|iVT(Z),,„l)|=i^E(^). 

Therefore, K^(z) controls the growth rate of the expectation of the number of 
particles which end up near z on the rescaled space. The next theorem shows that 
the supremum defining K^{z) corresponds to a unique optimal path hz] optimising 
over z then gives the total expected population growth. 

Theorem 4: 

For z G M, the optimisation problem 

K^{z) = K{hz,l) 

has a solution hz G C^[0, 1] which is unique for z ^ amongst all Hi functions 
ending at z. For z > 0, the solution hz is positive and satisfies for all s G [0, 1] 

h'^is) + m^phz{sy^^ = 0, hz{0) = 0, hz{l) = z. 



Furthermore there exists a unique % > such that 

Ke ■■= i^E(^E) = SUpi^E(^;) = sup K{f, 1) 
z /eC[0,i] 

Then the expected total population size satisfies 

lim^logE|iV(T)|=KE, 



T-5-00 rpj^ 



where one finds 



h',^{l) = 0, % = / ^\ , = (2m/3)-p 



/c 



1 dx 

v^T^^a^ 



2-p 






2 
2-p 



and 



Ke = —^m(3z^. 
2 +p 

Remark: 

For z < 0, one has hz{s) = —h-z{s). For z = and p > 0, there are two symmetrical 
optimal paths, one positive and one negative. For z = and p = the optimal path 
is unique and equal to Hq = 0. 

2.2 Almost sure growth along paths 

Let us now focus on the problem of giving an almost sure result for the actual 
number of particles that have a rescaled path lying in some set D. 
We let 

eoif) := inf {t G [0, 1] : K{f, t) < 0} G [0, 1) U {oo}. 

We think of Oq as the extinction time along /, the time at which the number of 
particles following / hits zero: if i > 9o{f), basically at large times no particle has 
a path that looks like / up to time t. On the other hand, if i < 9o{f), the number 
of particles with a rescaled path looking like / up to time t grows like the expected 
number of particles following that path. This is made precise in Theorem [5] below. 

Theorem 5: 

For any closed set D C C[0, 1] and t G [0, 1], 

limsup — 2+F ^°S I-^t(-D, f)| < sup{K{f,t) : f G D,6Q{f) > t} almost surely, 

and for any open set A C C[0, 1] and t G [0, 1], 

liminf ^2:p^log|A^T(-4,t)| > sup{K{f,t) : f G A,0o{f) > t} almost surely. 



Moreover, if one defines 

iCa.s.(^):=sup{K(/,l),/GC[0,l],/(l) = z,^o(/) = oo}, (15) 

we obtain the following corollary: 

Corollary 6: 

For zGR, 

lini lim — 5+^ log |-^t(^2 e, 1)1 = ^as (z) almost surely. 
e^or^oo T-ij^f 

Therefore, K^^s.{z) controls the growth rate of the almost sure number of particles 
which end up near z on the rescaled space. The next theorem shows that the 
supremum defining /Ca.sX-^) corresponds to a unique optimal path Qz (that, therefore, 
most particles ending up near z must have followed); optimising over z then yields 
the almost sure total population size growth. Let 

9 \ ^- / . 1 



r{s) := \-^{2 - pY \ and z := r{l) = \-^{2 - p)' \ . (16) 

Observe that by Theorem 111 for all s G [0, 1], Rst/T'^' '^~^' — ^ ?'(s) almost surely as 
T —7- 00. This means that r{s) describes the boundary of the limiting shape of the 
trace of the rescaled BBM and z is the rescaled position of the right-most particle 
at time 1. 

Theorem 7: 

For each z G [—z,z], one has 

K^.s.{z) = snp{K{f, 1) : / G C[0, 1], /(I) = z, \f{s)\ < r{s) Vs G [0, 1]}. (17) 

Moreover, the optimisation problem 

K^.,.{z) = K{g,,l) 

has a solution g^ which is unique for z 7^ amongst all Hi functions ending at z 
such that 9q{-) = cx). For \z\ > z one has ifa.sX-^) — ~oo, which means that no 
function of Hi with 9o{-) = 00 reaches z. 

The solution gz for < z < z is characterised as follows: there exists a unique 
Sz G [0, 1] such that 

(i) for all s G [0, Sz], gz{s) = r{s) 
(a) for all s G (sj;, 1], gz is twice continuously differentiable and 

g'^is) + m^pgzisy-^ = 0, gz{l) = z; (18) 

(Hi) gz is differentiable at Sz- 



Furthermore there exists a unique Za.s. > such that 

^a.s. := i^a.s.(^a.s.) = SUp K,,,,{z) = SUp {K(/, 1), / G C[0,l],eo{f) = Oo}. 

z 

Then the almost sure total population size satisfies 

lim — ^q:^ log I iV(T) I = Ka.s. almost surely, 

where one finds 

2 
2-p 



51,(1) =0, Za.s. 



y/2^ 



9 Iv / (JlJu 



+ 



^~P J2-^/P \/l - XP 



and 



2 — p 



Remark: 

It is easy to see that for z > we need only consider positive functions since for a 
general g,K{\g\,t) = K{g,t) for all t S [0, 1], and it is not hard to see that for all 
t, z > 0, Qzit) > 0. For p = 0, one has s^ = Vz G [0, z) and the almost-sure and 
expectation paths coincide. When p > 0, the proofs will make clear that Sz > 
Vz G [0, z]. In particular this means (still when p > 0) that the majority of particles 
found near the origin have in fact followed either the left-most or right-most path 
for some proportion of their history and travelled a long way out before increasing 
in number whilst heading back away from the frontier. 

An easy consequence of the theorem is that r(s) describes not only the limiting 
trace of the BBM, but also the actual rescaled trajectory of the rightmost particle at 
time T. It is the trajectory on the onset of extinction, the one for which K(r, t) = 
for all t or, equivalently, for which 

-r'{sf = m(3r{sf, with r(0) = 0, (19) 



as can be directly checked from (16). It is interesting to observe that r thus satisfies 
r"{s) = m/3pr{sy~^ . The solution g^ thus satisfies the same second-order differential 
equation on [0, s^) and (s^, 1] up to a sign difference on the second term. 

2.3 Probability of presence 

If / is such that 9o{f) < 1, Theorem Is] suggests that as T becomes large the number 
of particles whose rescaled paths have stayed close to / is almost certainly 0. The 
following large deviations result shows how the probability of presence of a particle 
close to / up to time t decreases as t goes from to 1. 



10 



Theorem 8: 

For any closed set D C C[0, 1] and t G [0, 1], 

limsup -^^ logF{NT{D,t) / 0) < sup 
r-s>oo 7^2:^ /eD 

and for any open set A C C[0, 1] and t G [0, 1], 



lim inf 

T-s>oo 



1 



2+p 



logF{NTiA,t) /0) > sup 



infi^(/,.) 



infK(/,s) 

s<t 



The case p = was proved by Lee [25] and again by Hardy and Harris 



2.4 Relationship to differential equations 

In this section we try to show how several of our results can actually be guessed 
from heuristic manipulations of partial differential equations. Although the whole 
discussion is informal it leads us to a theorem which gives an alternative description 
of Ka.s. and Ke- 

The expected density p{x, T) of points at position x and time T in the BBM we 
are studying can be written as the solution of the partial differential equation 

dp 1 d'^p 
df ~ 29x2 
Corollary [3] suggests that for large T, 



+ mj3\x\^^ p. 



(20) 



iogp(x,r) 



2+p 

T^-pK^{z) with z 



2 

T2-P 



(21) 



If we then plug (21) into (20) we get a differential equation: for T large, neglecting 
a term of order T^'^+p^ ^^ K^{z) , we get 



K^iz) + 



2z 



P 



2+p 
2-p 



K^iz) + 



2z^ 



{2-pf 



m(3\z\ 



(22) 



It is not obvious at first that the i^E(-z) defined in (14) is indeed a solution of (22) 



but we will show that this is the case in Theorem [9] below. The differential equation, 
however, is not enough to fully obtain K-e{z) as there is no obvious initial condition. 

A natural question is now: is there a differential equation of which i<'a,s.(-2) is 
solution? -fCa.s. {z) describes the growth rate of the almost sure number of particles at 
rescaled position z and it is different from the expected growth rate K¥,{z) because 
of some extremely rare events (on which particles go far away and reproduce a lot) 
which contribute to K¥.{z) and not to iCa.s. (-z) in the T — )• 00 limit. 

With this in mind, we now consider the inhomogeneous BBM with killing, where 
we remove any particles that ever cross the two space-time boundaries (s, ibx(s))s>o, 
for some given function x(-). The expected density p{x, T) of particles at position x 
at time T in this BBM with killing is then a solution of 



25x2 



+ rn,/3|x|^p, 



^p(±x(r),r)=o. 



(23) 



11 



If the absorbing boundary x(-) is taken to be the typical trajectory of the right- 
most particle, r(-), (see Theorem [I] and equation (16)), this in effect kills all those 
rare difficult paths and one might hope heuristically that p{x, T) describes, in some 
sense, the almost sure density of points in the original problem without the absorbing 
boundaries: 

logp(a;,T)~T^i^a.s.(-z) with z = ^^. (24) 

In fact, this heuristic does turn out to be the case. Indeed, the a.s. growth rate, 
I^a..a.{z)i is given by maximising over paths that end at z and do not undergo extinc- 
tion at any point, as in (15) and Corollary l6J However, ifa.s.(-^) can also be expressed 
by maximising over paths that end at z and never go beyond the right/left-most 
paths, as in (17) of Theorem ItI This equivalent representation is exactly the same 
as would be obtained for the expected growth rate in the BBM with killing at ibr(s), 
hence i^a.sX-^) is indeed given by equation ( p4| ) where p satisfies PDE (23). 

Further, we might even hope to determine r(-) in a self-consistent way: if, in (23), 



x{T) is significantly smaller than the almost sure position of the right-most particle, 
then many particles will gather close to the line and we can expect p{x{T) — 1, T) to 
be large. On the other hand, if x{T) is chosen significantly larger than the almost 
sure position of the right-most particle, then very few particles should come close to 
the boundary and p{x{T) — 1, T) should be small. Only for x{T) close to the almost 
sure position of the right-most particle can we expect p{x(T) — 1,T) to be of order 



one. We therefore reformulate (23) into 



( dp _1 d'^p 
df ~ 2 9x2 
p{±x{T),T)=0, 

dp 



+ m(3\x\^p, 



(25) 



dx 



{±x{T),T) 



±1, 



where we solve now simultaneously for the two unknowns p and x. 



We now plug (24) into (25) and, as the equations in the bulk for p and p are 
the same, we obtain the same equation (22) for i^a.s.(-2) as for Kf:^{z), albeit with 
different boundary conditions. From (25 ) we see that p{x, T) = for x > x(T). Thus 
(24) means that Ki^s(z) is not defined for z above some threshold value. Therefore 



we look for a solution of (22) only defined up to a finite value of z, which can only 



be the rescaled almost sure position z of the right-most particle, as defined in (16). 



The only way for this to happen is for the right-hand side of (22 ) to vanish at z 



z: 



2+p 
2-p 



K^.s.{z) + 



2z' 



2-py 



ml3\z\- 



0. 



(26) 



(When the right-hand side reaches 0, one can check that the second derivative di- 
verges and the solution cannot be continued beyond that point.) Furthermore, 
one must have i('a.s.(^) = from the second boundary condition in (25); indeed 



Ka.,s.{z) > would correspond to having increasingly many particles at the bound- 
ary while i^a.s.(^) < would mean that the number of particles next to the boundary 



goes to zero. We thus recover the expression for z given in (16) as a solution of (26). 
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Prom the descriptions of K^[z) and i^a.s.(^) given in Theorems H and ^ one can 
show that, indeed, 

Theorem 9: 

K^{z) and Ks_^s{z) as defined in Theorems^ and^ are, for z > 0, two solutions of 
the following differential equation: 



K'{z) 



2z / 2 + p^^, , 4z2 



P 



P 



{2-pf 



2m(3zP. 



(27) 



(Note that (27) is not imphed by (22) as we could have put a minus sign in front 
of the square root. We will show in the proof section that K'{0'^) > which justifies 
the choice of the plus sign. For 2: < 0, by parity of K{z), the other sign must be 
chosen.) 

This approach from partial differential equations can be extended to the case 
p = 2. Considering only the particles that do not go further away than the almost 
sure position of the right-most, we start from (25) with p = 2. The scaling function 

(28) 



(24) obviously does not work for p = 2, but if one plugs 

with z = xe~ 



log p{x,T) 



e'^^L{z) 



into (25) one gets that 
1 



1 



-yL'{z) + AzY = 2AL{z) + -A^z^ - m(3z'^ 



(29) 



where a term of order e~'^ L" (z) has been neglected. (29) has exactly the same 
structure as ( 22 ) except that yl is a priori an unknown quantity. However, requesting 
as for p < 2 that the right-hand side vanishes at z = z and that L{z) = implies 
that 

A = y^2mP (30) 



so that the equation reads 



l{L'{z) + ^/2^z)^ = 2^/2^ Liz) 



(31) 



or, taking the square root for z > 0, 

L'{z) = -^J2mj3z + 2\l ^j2ml3L{z). (32) 

(We put a plus sign in front of the square root by analogy with the p < 2 case.) 

Note already that ( 30 ) with ( 28 ) allows one to recover the results from Berestycki 
et al. [5j: 

For p = 2 and m = 1, lim — log log |A^(T)| = 2y 2/3 almost surely, (33) 

T->oo 1 

and from Harris and Harris |16j : 

For p = 2 and m = 1, lim — log-R^ = \/2/3 almost surely. (34) 

T— >-oo 1 



13 



(Although the two papers [5j and [16] only concerned the binary branching (m 



1) case, their results (33) and (34) could easily be extended to the more general 
branching process of the present paper, with an arbitrary value of m.) 



Going further, (32) can be solved by making the change of variable L{z) 
^J2m|3 z^(l){z)'^ with cl){z) > 0. One gets for ^; > 



2z(l)'^{z) + 2z^(l){z)(l)'{z) 
For z > the variables can be separated 

2cj,{z)cP'{z) 



2(/)2(z) - 2(j){z) + 1 



-z + 2z(i 

1 

z 



(35) 



(36) 



and, after integration of both sides and simplification, one gets an implicit form for 
Liz): 



L{z) = y^2mf3z'^^^ 
exp [arctan(l 



C- 



m 



with 6 > 



^^^^^ ^^ (37) 

^202 _ 20 + 1 """ "^ " " 

where C is an integration constant. By taking (/> — t- oo one obtains L{0). The 
rescaled position z of the right-most particle is when L{z) = or (j) = 0. From 
(31), the optimal position z where L is maximal is such that L(z) = ^J2mfi z^ /A or 
(j) = 1/2. This leads to 



L(0) 



(J2 

^/2^—e 



Ce5, z = CV2, L{z) = y/2^ 



C 



(38) 



Remark: 

We have no theory for the value of the integration constant C; in fact we be- 
lieve that C is realisation-dependent. This method allows us however to make 
some conjectures on the values of several ratios such as zjz =(the position of 
the right-most) /(the position where the density of particles is the highest), or 
L{z)jz^ = log I A'r I /(position of the right-most )2. We have no demonstration for 
these conjectures. We simply observe that, for instance, the value of the ratio zjz 
computed in the p <2 case (see equation ( 16 ) and Theorem [7]) converges as p — )• 2 
to the ratio predicted by ( 38 ) . 



2.5 Explicit calculations for p = 1 

It is interesting to note that for p = 1 (as well as for the easier case p = 0) the 
equations given by Theorems |4] and [7] can be solved explicitly. 

• The rightmost particle satisfies Rt ~ j"^/^^'^ ^■^■ 

• The optimal path for expected growth with end-point z > is given by 



with growth rate 



1 1 

hz{s) = --m(5s^ + {z+ -ml3)s 



Kdz) = ^"1^/3^ + 2"^/^^ - 2^^- 
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The optimal end-point for expected growth is % = j'ti/S, giving the optimal 
path 

/i%(s) = --m/3s'^ + ml3s 
and total expected growth rate Ke = ?n.^/3^/6. 
r{s) = ^mjis^ and z = ^rn,/3. 
The optimal path for almost sure growth with end-point z G [0, 2] is given by 



9z{s) 



-'^s'^ + 2m(3szs-m(3sl ifsG(s^,l] 



where s^ = 1 — ^/i ^. The corresponding growth rate is 



2 m/3 

• The optimal end-point for almost sure growth is ia.s. = 41 giving the optimal 
path 

,,_r^.2 if^e[o,i] 

5£a.e.l5j-| _r^^2^^^^_m^ if5G(i,l] 

and total almost sure growth rate -fCa.s. = m^/3^/12. 

We note in particular that there is positive expected growth for all \z\ < m/3(l/2 + 
l/\/3) despite almost sure extinction for all \z\ > m/3/2. Only for p = do the 
almost sure and expected growth match. 

2.6 Proof strategy 

As already pointed out in the introduction, the results in expectation (Theorems 
P^ and El) can be derived in a rather straightforward fashion from Shilder's large 
deviation theorem and the use of the so-called many-to-one principle (see SectionH|. 
Some fairly standard large deviations techniques (using for example the exponential 
tightness of Brownian motion) can then be used to deduce the large deviations 
behaviour of the system seen in Theorem Is] (see Section IgI) . 

For the almost sure growth along paths, however, we need something stronger. 
Using the many-to-one lemma we construct in Section [3] processes that are non- 
negative martingales which count the numbers of particles whose paths lie in certain 
sets. We then use the fact that these martingales have almost surely finite limits 
to obtain an almost sure upper bound on the number of particles whose rescaled 
paths remain close to /. It is also quite usual, at least within the world of branching 
processes, that if an additive martingale — like ours — is uniformly integrable, then 
it has a strictly positive limit, giving us our almost sure lower bound. Again this 
is true in our case, although showing it is highly non-trivial — a large part of the 
work for this article is spent in proving this lower bound. 
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Figure 1: For p = 1 and mf3 = 1, the thick (red) line in the left-hand graph shows 
r(s), the (rescaled) path taken by the right-most particle. The thinner lines show, 
for various values of the endpoint z, the optimal paths Qzis) for almost sure growth 
(plain blue lines) and hz{s) for expected growth (dashed green lines). In the right- 
hand graph we show the profiles i^a.s. (-z) and Ke{z) of the number of particles alive 
at each position z for almost sure growth (plain blue) and expected growth (dashed 
green) . 

We then set about proving the results concerning how many particles follow 
particular paths. In Section [5] we prove Theorem [5j applying many of the results 
obtained in the previous two sections. 

We move on in Section [7] to derive the optimal paths seen in Theorem [7] and 
Theorem Hi and study these paths further (in particular proving Theorem ^ in 
Section [H 

3 A family of spine martingales 

3.1 The spine setup 

A key idea in our proofs will be the use of certain additive martingales. These 
martingales can be used to define changes of measure under which one particle 
behaves differently than under the law P of our branching particle system. The 
tools introduced in this way are extremely useful, and should be fairly intuitive. As 
they are now well-embedded in the branching process literature, we will leave out 
several proofs in this section, and refer the interested reader to Hardy and Harris' 
very general formulation in ^15j. 
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We first embellish our probability space by keeping track of some extra infor- 
mation about one particular infinite line of descent or spine. This line of descent is 
defined as follows: our one initial particle is part of the spine; when this particle dies, 
we choose one of its children uniformly at random to become part of the spine. We 
continue in this manner: when the spine particle dies, we choose one of its children 
uniformly at random to become part of the spine. In this way at any time t > 
we have exactly one particle in N{t) that is part of the spine. We refer to both this 
particle and its position with the label ^t] this is a slight abuse of notation, but it 
should always be clear from the context which meaning is intended. The spatial 
motion of the spine, {Ct)t>o, is a standard Brownian motion. 

The resulting probability measure we denote by P, and we find need for four 
different filtrations to encode differing amounts of this new information: 

• J-'t contains all the information about the original system up to time t. How- 
ever, it does not know which particle is the spine at any point. Thus it is 
simply the natural filtration of the branching Brownian motion. 

• J^t contains all the information about both the BBM and the spine up to time 
t. 

• Qt contains all the information about the spine up to time t, including the 
birth times of other particles along its path and how many children are born 
at each of these times; it does not know anything about the rest of the tree. 

• Qt contains just the spatial information about the spine up to time t; it does 
not know anything about the rest of the tree. 

We note that J^t C J-t and Qt d Qt '^ ^t, and also that P is an extension of P in that 

Lemma 10 (Many-to-one lemma): 

U 9it) ^s J-t-measurahle it can he written 

ait) = Yl 9uit)l{^,=u} 



where each gu{t) is J-t-measurahle, and then 



E 



ueN{t) 



9u{t) 



]E[e"^/ol«»l^ds^(t)]. 



This lemma is extremely useful as it allows us to reduce questions about the entire 
population down to calculations involving just one standard Brownian motion — the 
spine. A proof of a more general version of this lemma may be found in [15| . 
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3.2 Martingales and changes of measure 

For p £ [0, 2), / G C[0, 1], G [0, 1] and e > 0, let 

2 

q:= G l,oo) 

2 — p 

and define 

NT{f,e,9) := {u G iV(^T) : |X„(t) - T'^fit/T)\ < eT^ Vt G [0,er]} 



so that NT{f,e,e) = NT{B{f,e),9), see ([13J), where 

B{f,e):={geC[0,l]:\\f-g\\oo<e}. 

We look for martingales associated with these sets. For convenience, in this section 
we use the shorthand 

NT{t):=NT(,f,e,t/T) 

and 

CT(x,t):=cos(^(x-TV(t/T))). 

The following Lemma is adapted from Lemma 6 in |I8j . 

Lemma 11: 

If f e C"^ [0, 1] then the process 

is a Qt-local martingale under P. 

Proof. Since the motion of the spine is simply a standard Brownian motion under 
P, this is easily checked by applying Ito's formula (the sufficient conditions of, for 
example, Lawler ^24j tell us that if / G C^[0,1] then Vt is sufficiently smooth for 
Ito's formula to hold). See |18j . Lemmas 5 and 6, for the calculations. D 

By stopping the process (^^(i),* G [0, T]) at the first exit time of the Brownian 
motion from the tube {(x,i) : \T'^f{t/T) — a;| < eT'^}, we obtain also that 

C,T{t) := VT{t)li^\Tif{s/T)-£,s\<eTi Vs<t}, * ^ [0, T] 

is a non- negative ^^-local martingale, and since its size is then clearly constrained 
it must in fact be a ^(-martingale. As in [15], we may build from C,t a collection of 
J^i-martingales C,T given by 

Qrit) := n(l + ^.)e-"^^«l«=l'"^Cr(t), tG [0,T]. 

When we project CTit) back onto J-t we get a new set of mean-one J^i-martingales 
{ZT{t),t > 0). These processes Zt are the main objects of interest in this section, 
and can be expressed for t £ [0, T] as the sum 

ZTit)= Yl C^^\t)e-"'^fo\^-i-^r'^^= Y^ yH(^)g-m/3/„^|X4.)|-ds 
uGNrit) ueNrit) 
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where Ct (^) ^^"^ V^" (t) are simply Crit) and Vrit) with the path Xu{s) replacing 
C{s), i.e. 

We now proceed to show that the martingales Zt are close to 

Ehu. close to /le^"-^ ^0 [|/'(^/T)-™m(VT),] d. (33) 



and that they have the properties that we discussed in Section 2^ — specifically, 
we aim to show that for certain / the martingales Zj- are uniformly integrable and 
thus cannot be too small. This is the key step to counting particles whose rescaled 
paths stay close to /. 

We define new measures, Qt, via 

QtI^, := Cr(t)lP|^, 
for t E [0, r] — and note that 

QtIj-, = ^T(t)lPb, and Qrh = CrirnGf 

Lemma 12: 

Under Qt the spine (i^i,t G [0;^]) nioves as a Brownian motion with drift 

T'^-'fit/T) - ^ tan (^(x - T'^fit/T))) 
when at position x at time t; in particular, 

16 - T'if{t/T)\ < eT'i yt < T. 

Each particle u in the spine dies at an accelerated rate {m + l)f3\x\P when in position 
X, to be replaced by a random number A^ + 1 of offspring where A^ is taken from the 
size-biased distribution relative to A, given by Qt(^u = k) = [m + l)~^{k + l)P{A = 
k) (note that this distribution does not depend on T). All non-spine particles, once 
born, behave exactly as they would under F: they move like independent standard 
Brownian motions, die at the normal rate f3\x\P, and give birth to a number of 
particles that is distributed like 1 + A. 

Proof. A proof of this result can be found in ^15^ . We will not use the precise drift 
of the spine except for the fact that it remains within the tube: to see this note that 
since the event is ^T-measurable, 

Qr(3t < T:\Ct- T^fit/T)] > eT'^) = nQT{T)l{^t<T:\i,-T^f{t/T)\>eT^)] = 

by the definition of C,t{T). □ 

Another important tool is the spine decomposition, which will allow us to bound 
the growth of the martingales Zt via one-particle calculations in a more delicate 
way than is possible via the many-to-one lemma. A proof of a more general version 
of the spine decomposition may be found in |15j . 
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Theorem 13 (Spine decomposition): 

Q^T-o-linost surely, 

qT[ZT{t)\QT\ = Yl ^«VT(5„)e-™'^io'" l««l^ds ^ y^^^y-n^^J^li^lPds 

where {u < S,t} is the set of ancestors of the spine particle at time t, and Su denotes 
the time at which particle u died and split into 1 + A^ new particles. 

As we have aheady mentioned, the main aim of introducing these martingales is 
to give us a lower bound on the number of particles in Nxit). To do this we must 
bound the size of each of the terms in the sum. The following lemma is a simple 
bound for the Girsanov part of the martingale, using integration by parts. 

Lemma 14: 

/// G C^[0, 1] and /(O) = then for any u G Nxit), almost surely under both IP and 
Qt we have 

r«-i f f'{s/T)dXu{s)-T^'^~^ f f'{s/Tfdt 
Jo Jo 

<2er29-2 f \f'\s/T)\ds + eT^''-^\f'{0)\. 
Jo 

Proof. ^From the integration by parts formula for Ito calculus (since for any particle 
u £ N(t), {Xu{s),0 < s < f) is a Brownian motion under P) we know that for any 
ge C2[0, 1], under P, 

g'{t)Xu{t)= [ g"{s)Xu{s)ds+ I g'{s)dXu{s). 
Jo Jo 

^From ordinary integration by parts, if ^(0) = 0, 



t rt 



g'{syds = g'{t)g{t)- / g{s)g"{s)ds. 
Jo 

Now set g{t) = r«/(t/T) for t G [0,T]. We note that if u G Nrit) then |X„(s) 
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g{s)\ < eT'i for all s<t. Thus 

T«-i f f'{s/T)dXu{s)-T^'^-^ f fis/Tfds 
Jo Jo 



< 



g'{s)dXu{s) - [ g'isfds 
Jo 

g'{t){Xu{t)-g{t))- [ g"is){X^{s)-gis))ds 



< \g'{t) - g'{0)\ X \X^{t) - g{t)\ + \g'{0)\ x |X„(t) - g{t)\ 

+ [ \g"{s)\x\Xu{s)-g{s)\ds 
Jo 

<2eTi I \g"{s)\ds + £T'i\g'{0)\ 
Jo 

= 2er25-2 f |/"(s/r)|ds + er29-V(0)| 
Jo 

almost surely under P and, since Qt ^ IP (on J-t), almost surely under Qt- D 

The next lemma continues along the same theme, controlling the terms in Zt 
so that eventually we will be able to give a lower bound on Zt and then use this to 



give a lower bound on Nt by showing that Zt looks something like (39) 



Lemma 15: 

For any u G NT{t), 



T29-2 inf f \g{s/T)\Pds< [ \Xu{s)\Pds<T^''-^ sup [ \g{s/T)\Pds. 

9<^B{f,e) Jo Jo geB(f,e) Jo 

Proof. This follows immediately from the fact that if n G NT{t) then (by definition) 
there exists g S B{f,£) such that Xu{s) = T'^g{s/T) for all s <t. D 

We are now ready to start putting together the bounds that we have given. 



Combining Lemmas 14 and 15 with the spine decomposition we obtain the following. 



Lemma 16: 

/// G C^[0,1], /(O) = 0, /'(O) = and m(3 J^ \f{s)\Pds > ^J^* /'{s^ds for all 
(p G (0, 0], then for small enough e > and any T > and t < 6T, there exists ry > 
such that 

QT[ZTit)\QT] < V ^„e^'/(S^'^''~')-''"^'^i'o'" |6l^ds ^g7rV(8e2r2«-i)-^^m/3/o*|6|^d. 
w<6 



-almost surely. 
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Proof. Recall that under Qj' the spine is in Nxit) for all t < T. Thus by Lemmas 



14 and 15, since /'(O) = 0, for any 77 > 

-m(3 Tie^rds + r"-! ff'{s/T)di,-l-T^'^-^ f f'is/Tfds 
Jo Jo ^ Jo 

<-r]mf3 I |6rds-(l-7?)m/3r2'?-2 inf / \g{s/T)\Pds 

Jo <J&B{f,e) Jo 

t 



^^r^2g^2 f f(^s/Tfds + 2eT^'^-^ [ |/"(s/T)|d 
2 Jo Jo 



for all t < T. Then, since m/3 f^ \f{s)\Pds > \ f^ f'{sfds for all </> G (0, 9], we may 
choose e > and 77 > small enough such that 



;i - 7?)m/3r2'?-2 inf / \g{s/T)\Pds 

9<^B{f,e) Jn 



^1^2g-2 f /(s/T)2ds + 2er2'?-2 [ |/"(s/T)|ds < 
2 Jo Jo 

for all t G [0, OT]. Plugging this into the spine decomposition, we get 

We are in a position now to complete one of our two initial aims, which was to 
show that the martingales Zt are uniformly integrable. This will be used in Section 



5.1 to show that Zt cannot be too small. 



Proposition 17: 

Take f G C^[0,1] and 9 G [0,1]. // /(O) = 0, /'(O) = 0, and mP J^ \f{s)\Pds > 

2 /o /'(•5)^ds for all (j) G (0,0], then for small enough e > the set {Zxit) : T > 
1, t < 9T} is uniformly integrable under P. 

Proof. Fix 5 > 0. We first claim that there exists K such that 

sup QT{QT[ZT{t)\GT] >K)< 6/2. 
T>1 

t<eT 

To see this, take an auxiliary probability space with probability measure Q, and on 
this space consider a sequence Ai,A2,... of random variables with the same (size- 
biased) distribution as A under <Qt (there is no dependence on T) and a sequence 
ei , 62 , . . . of random variables that are exponentially distributed with parameter 
f3{m + 1); then set 5„ = ei + • • • + e„ (so that the random variable Sn has the same 
distribution as Jq " l^^l^ds, where Su is the time of the nth splitting event along the 
spine under Qr)- By Lemma 16 we have (since 2*7 — 1 > 1) 

(00 
Y^A.e^'/''''-'^^^ +e^'/^'' >k] . 
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Hence our claim holds if the random variable 






;e 



-r)Sj 



can be shown to be Q-almost surely finite. Now for any 7 G (0, 1), 






g-r?5„ = 00 < Q{Ane-''^" > 7" infinitely often) 



< Q { > log 7 H infinitely often 1 . 

n n 



By the strong law of large numbers, Sn/n — >• 1/P{m + 1) almost surely under Q; so 
if 7 G (exp(— r///3(rn- + 1)), 1) then the quantity above is no larger than 

/ log^n ^ „\ 

Q I lim sup > (J I . 

\ n— ^-oo IT- / 

But this quantity is zero by Borel-Cantelli: for any T, 

n ^ ^ n 

< / QT(logA> ex)dx 
Jo 

'log A' 

e 

which is finite for any e > since (by direct calculation from the distribution of A 
under Qr given in Lemma 12) Qr[log^] = P[Alog A] < 00. Thus our claim holds. 
Now choose M > such that 1/M < 5/2] then for K chosen as above, and any 
T>l,t<dT, 

QriZTit) > MK) < QTiZrit) > MK, QT[ZT{t)\GT] < K) 

+ Qr(QT[^T(i)|aT] >K) 
'ZtH). 



< 



MK '^{^T[ZT{t)\gT]<K} 

QT[ZT{t)\GT] 



+ 6/2 



-1, 



lT[ZT(t)\gT]<K} 



MK 

< 1/M + 6/2 < 6. 

Thus, setting K' = MK, for any T > 1, t < 6T, 

F[ZTit)l{Zr(t)>K'}] = ^HZrit) > K') < 6. 
Since 6 > was arbitrary, the proof is complete. 



+ 6/2 



n 
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Finally we show that Zq^ is close to (39) 



Lemma 18: 

For any 5 > 0, if f & C^[0, 1], /(O) = and e is small enough then 



ZT{eT)<\NT{f,e,e)\e^^' 



nO_^pT^,-l jO |^(,)|Pd,+ 1^2,-1 J9 ^,(^)2d^+5^2,-l 



Proof. Simply plugging the results of Lemmas 14 and 15 into the definition of Zt{OT) 
gives the desired inequality. D 

We note here that, in fact, a similar bound can easily be given in the opposite 
direction, so that \NT{f,£/2,6)\ is dominated by Z^{9T) multiplied by some deter- 
ministic function of T. We will not need this bound, but it is interesting to note 
that the study of the martingales Zt is in a sense equivalent to the study of the 
number of particles Nt- 



4 Proof of Theorem [2] 

We first rule out the possibility of any particles following unusual paths, which allows 
us to restrict our attention to compact sets, and hence small balls about sensible 
paths. In this section we go back to our original notation convention in which small 
letters such as s and t as well as 9 are used for scaled time parameters varying in 
[0, 1] and capital letters are reserved for the non-scaled time. 

Lemma 19: 

Fixde [0,1]. For N en, let 

Fn := |/G C[0,1] ■.3n>N,s,t£ [0,9] with \t - s\ < ^, |/(i) - /(s)| > ^ 

For any r] > we may choose N gN such that for all large T 

F{NTiFN,9) / 0) < E[Nt{Fn,9)] < eM'VT^'^'^) 
and 



limsup -2— Y log {NxiFN, 9] 

almost surely. 

Proof Fix T > 5 > 0; then for any U G [S, T], 

{^U e Nu{Fn, 9)} = hn >N,s,te [0, 9]:\t-s\< 

C \3n>N, s,t£ [0,9] : |t-s| < 



-oo 



n^ 



^sU - S,tu 



n^ 



6t - CtT 



Si 



> 



> 
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Since the right-hand side does not depend on i, we deduce that 

{3Ue[S,T]:^u(^Nu{FN,9)} 



1 



C <3n>N, s,tG [0, 6*] : |t - s| < -^, 



n^ 



S,sT - itT 



Si 



> 



Now, for s G [0,9], define 7r(n, s) := [2n^sJ/(2n^). Suppose we have a continuous 
function / such that sup^grgei \f{s) — f{7r{n,s))\ < l/{4:^/n). If s,t £ [0,9] satisfy 
|t — s| < 1/n^, then 

1/(0-/(^)1 
< \fit) - finin, t))\ + \f{s) - fiirin, s))\ + |/(7r(n, s)) - finin, t))\ 



1 1 2 

< — = + -^ + 



1 



6t - ■^7r(n,s)T 



Si 



> 



1 



40iJ ■ 



4-^/n 4-y/n 4-y/n ^/n 
Thus 

Standard properties of Brownian motion now give us that 

F{3U G [S,T] : Cf/ G Nu{Fn,9)) <F{3n>N,s<9: |6t - e.(n,s)Tl > SyA^n) 

< 



V 2n¥ [ SI 

n>N \*e[0'l 



sup |6t| > 5V4\/^ ) 

/(2n2)] y 



16T 



n>7V ^ ^ 

Taking S = j and T = j + 1, we note that for large N, 






n>Af 



n>N 
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Now, for any M > 0, 



sup \Nt{Fn, 

.T€[j,j+1] 



> 1 <E 



sup \Nt{Fn, 
Te[j,j+i] 



< E 



< 



Z2 ^{^Te[j,j+i]:ueNT(FN,e)} 
ueN(j+i) 



E 



1; 



™/9/o I6l''ds| 



{3Te[j,j+mTeNT{FN,0)} 



{3Teb,j+l]:5TGiVt(Fjv,e)}J^{sups<^.+i|?sl<M(i+l)n 



+ E 



m/5j^^'\U^ds-^ 



{suPs<j+l l€sl>M{j+l)9} 



< e™^^''(^-+i)^^"V(3T G [j,j + 1] : Ct G iVT(F^,e)) 



fc>i 



{sups<j+i\i;T\e[kM{j+l)i,{k+l)MU+m} 



,m/3{j+l)29-i(fe+l)PAfPp 

< e-^^''(^-+i)^^"V(3T G [j,j + 1] : Ct G NriF^^e)) 

A:>] 



+ ^ gm/^U+i)- n'c+irM''p( g^p l^^l g [^^(j- + ^y^ (^ + i)^(^- + 1)9]) 



^ V27r(i + 1) 



Both of the terms in the right-hand side can be made exponentially small in j by 
choosing M, and then A^, sufficiently large (for the first, see our calculations earlier 
in the proof). This establishes the first part of the lemma, and by Borel-Cantelli we 
have that for large enough N 



lim sup sup I Nt (-F/V , 
j^oc Telj,j+i] 



> 1) = 



and since \Nt{Fj^,6)\ is integer-valued. 



1 



lim sup -2— J log \Nt{Fn, 



-oo 



almost surely. 

We now check that we can cover our sets in a suitable way. 
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Lemma 20: 

Let 



Co[0,1]:={/gC[0,1]:/(0) = 0}. 



For each N £ N, the set Co[0, 1] \ F/v is totally bounded under \\ • ||oo (that is, it may 
be covered by finitely many open balls of arbitrarily small radius). 
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Proof. Given e > 0, choose n such that n > N V 1/e^. For any / G Co[0, 1] \ F/v, if 
|ti — s| < 1/n^ then \f{u) — f{s)\ < l/\/n < £■ Thus Co[0, 9] \ F/v is equicontinuous 
(and, since each function must start from 0, uniformly bounded) and we may apply 
the Arzela-Ascoli theorem to say that Co[0, 1] \ F^ is relatively compact, which is 
equivalent to totally bounded since (C[0, 1], || • ||oo) is a complete metric space. D 

Lemma 21: 

ForDc Co [0,1] define 

D' := {/ G C7o[0, l]:3g£D with \\g - /||oo < e}. 

For any 5 > there exists e > such that 

snpK{f,e)<svLpK{f,e) + 6. 
fGD'^ feD 

Proof. For each 0, f '-^ L f'{s)'^ds is a lower semicontinuous function on Co[0,^]: 
we refer to Section 5.2 of [H] but it is possible to give a direct proof. Thus / i— >• 
m/3 Jq \f{s)\Pds — 5/0 /'(s)^ds is upper semicontinuous. Now, by Jensen's inequal- 
ity, for any / S Co[0, 9] and any s,t £ [0,9], s < t, 

" ' 1 I'., >j \^ (nt)-j{,Y'' 



^Jj'(u)Mn>[^Jj'iu),u 



so that 



if{t)-f{s))^<{t-s) f fiufdu. (40) 

J s 



There exists t G [0, 9] such that \f{t)\P > I /^ \f{s)\Pds, so by (|40|) (taking s = 0) 

i /'W^dn > I f\uf du > ^ ^,^^^ ^ > (y^ |/(.)rd.) 
and hence 

{/ € Co[O,0] : m(3 [ \fis)\Pds -\ f f'^sfds > K} 
Jo ^ Jo 

/ f-e xP/2 . .e 

C {/ G Co[0, 0] : m/3 H /'(s)'ds J - - ^ /'(s)2ds > K} 

c{fGCo[0,9]: [ f'is)^ds<K'} 
Jo 



for some K' since p/2 < 1. But by (40), 



{/G Co [0,0]: /" /'(s)2ds<K'} 

JO 

C {/ G Co[O,0] : Vs,t G [0,0], |/(s) - f{t)\ < ^{t - s)K'} 
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and the Arzela-Ascoli theorem teUs us that this latter set is totally bounded. Thus 
the set 

{/ e Co[0, e]:ml3 [ \fis)\fds -\ f f'isfds > sup K(/, 9) + 5} 
Jo ^ Jo feD 

is totally bounded, but by upper-semicontinuity it is closed, and hence compact. 
Since it is disjoint from {/ G Co[0, 6] : 3g £ D with f{s) = g{s) Vs G [0, 9]}, which is 
closed, there is a positive distance between the two sets. D 

Before continuing, we need the following lemma. 

Lemma 22: 

For any x,y £M and p G [0, 2), 

\x + y\P < |x|P + |y|P + 2|x|P/2|y|P/2. 

This result is entirely elementary; see |27] for a proof. 

Proposition 23: 

If f ^ -F/V; then for A = B{f, e), we have 

limsup-2^1ogE[|ArT(A^)|] < sup K{g, 9) + RNie) 



and 



T^oo J^--' g^A 



hminf -^ logE[|iVT(A ^)|] > supi^(5,^) - i?7v(e) 



^^(^) ■= { n„.^/Af2+i , -\P/^ 



as T ^ OG, where 

ifp = 

2m/3(^ + e)^^^2e)P/2 + (2e)P tf p > 0; 

in particular R is a deterministic function of e such that for each N, -RAr(e) — )■ as 

Proof. ^From Schilder's theorem (Theorem 5.1 of [29]) we have 

limsup -^ logP(eT G Nt{A, 9)) < - inf J / f'{sfds 
and 

Thus, by the many-to-one lemma, 



1 , „r, ,, , T „s n . 1 - _ r ™a per I 



limsup-2— Y log E [I iVT(^, 6*) I J < lim sup -2— j log E 

T->oo -^ T— !>oo J 



<supm/3 / |(7(s)|^ds— inf - / g'{s) ds 
neA Jo g&AnHi 2 Jq 



geA 
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and similarly 

l^^rif -^logE[\NT{A,e)\] > inf m/? / \g{s)\Pds - inf \ [ g'{sfds. 



Note that since / ^ Fat, sup^gjo^^] 1/(^)1 < (^ + l)/v^ (split [0,6*] into N^ + 1 
intervals of equal width; then / changes by at most 1/yN on each interval). Now 
fix (5 > and choose h £ ACi Hi such that 

/ h'isfds< inf / g'isfds + 5. 

For any 5(1,5-2 e A, 

9i{sWds< [ {\g2{s)\+2eyds 
Jo 

p/2 



g2{s)\Pds + 2 [ \g2{s)\P/\2er/'ds+ [ {2eyds 
Jo Jo 

re / /V2 _L 1 \ P/2 

< J^ \g2{sWds + 2 i^—^ + ej (2efl^ + (2eYl^ 



and similarly 



5i(s)Rs> / |52(s)rds-i?jv(e). 

^0 



Thus 



inf m/3 / \g{s)\^ds- inf J / 5'(5)'d5 

,e -, /.(? 

>snpm^ \g{s)\Pds- mi - g'{sfds-RN{e 

geA Jo gdAnHi 2 Jq 

> sup \m^ [ \gis)\Pds-l [ g'isfds}-RNie) 
geAnHi I Jo '^ Jo 



and 



rd I [^ , 

supm/3 / \g{s)\^ds — inf - / g'{s) ds 
gdA Jo gf^AnHi 2 Jq 

pe 1 r-e 

<m/3 \h{s)\Pds-- h'{sfds + RN{e) + 6 
Jo 2 Jq 

< sup |m/3 / \g{s)\Pds-l I g'{sfds\+RN{e) + 5. 

g&AnHi I Jo ^ JO J 

Since 6 > was arbitrary, this gives the desired result. D 
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Proof of Theorem^ By Lemma 19 we may first choose A^ G N such that 



E[A^T(i^7V, ^)] < exp f r2«-i(sup ir(/, 0) - 1) J 

for all large T (just pick ry = supjgj^ K{f, 0) — 1 > 0). 

Let us now fix 5 > and (by Lemma [2l|) choose e > such that 



sup K{f, e) + RN{e) < sup K{f, 9) + 5. 



Then (by Lemma 20 ) choose n and /i , . . . , /„ such that 

n 

Ec(JB{fi,e)cE' 

i=l 

where for the lower bound we choose E to be an open set A, and for the upper 
bound we choose E' to be a closed set D. 

Let us first deal with the lower bound. For any f £ A and e > so that 
B{f, e) C A we can use Lemma [23| to obtain 

limmiT^-'^nogE[\NT{A,e)\] > liminf T^^^gj^g^^i^^^^^j^ ^^^ ^)|] 

> sup K{g,9)-RN{e) 

and since e is arbitrarily small we immediately see that 

liminfr^-2giQg]£[|^^(^^^)|] > K{f,e). 



Let us now turn to the upper bound. By Proposition 23 we have 



limsup-^-^logE[\NTiD, 



< sup K{f,e)+RN{e) 
feD^ 

< sup K(f, 9) + 6 

where the last inequality is due to our choice of N and e. Since (^ > is arbitrary, 
the result follows. D 

Proof of Corollary]^ Observe that if we take D = {/ e C[0, 1] : /(I) G [a, 6]} 
and A = {f £ ^[0,1] : /(I) € (a, ^)} {D is closed and A open) then necessarily 
sup j-^^ K{f,t) = sup jr^j^K{f,t). To see this observe that for any / such that, say, 
/(I) = a one can find a sequence /„—)•/ in ^ such that K(fn, 1) — )• K{f, 1) (just 
modify / near time 1 to finish ever so slightly above a). D 
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5 Proof of Theorem [5] 

5.1 The heuristic for the lower bound in Theorem [5] 

We want to show that NT{f,£,0) cannot be too small for large T. Recall that for 
/ e C[0, 1] and 6 e [0, 1], we defined 

K{f,e):=l rnl3f^\m\Pds-Uof'('?^^ ii f ^ H^ 
\ — oo otherwise. 

which is the growth rate in expectation. We are going to show that the almost 
sure behaviour is described by a rate function which differs by the presence of a 
truncation at the extinction time ^o- 

Step 1. Consider a small rescaled time rj. How many particles are in NT{f, e, r])? 
If r/ is much smaller than e, then (with high probability) no particle has had enough 
time to reach anywhere near the edge of the tube (approximately distance eT from 
the origin) before time rjT. Thus, with high probability, 

|iVT(/,e,r/)| = |iV(r?r)|. 

We can then give a very simple (and inaccurate!) estimate to show that for some 
u > 0, with high probability, 

\N{r]T)\ > vT. 

Step 2. Given their positions at time jyT, the particles in Nxif, £, rj) act indepen- 
dently. Each particle u in this set thus draws an independent branching Brownian 
motion. Let Nt{u, f, e, 9) be the set of descendants of u that are in Nxif, £, 0). How 
big is this set? Since ry is very small, u is close to the origin at time r]T. Thus we 
may hope, for a given 5 > 0, to find some 7 < 1 such that (for each u) 

F{\NT{u,f,e,B)\<eMK{f,e)T'''~'-6T'^~'))<^. 

Step 3. If NT{f,£,0) is to be small, then each of the sets Nt{u, f,e,9) for 
u £ NT{f,£,r]) must be small. Thus 



F{\NT{f,e,e)\<exp{K{f,e)T^''-'-6T^'^-')) < 



7^^, 



and we may apply Borel-Cantelli to deduce our result along lattice times (that is, 
times Tj, j > such that there exists r > with Tj — Tj_i = r for all j > 1). 

Step 4- We carry out a simple tube-reduction argument to move to continuous 
time. The idea here is that if the result were true on lattice times but not in 
continuous time, the number of particles in NT{f,£,0) must fall dramatically at 
infinitely many non-lattice times. We simply rule out this possibility using standard 
properties of Brownian motion. 

The most difficult part of the proof is step 2. However, the spine results of 
Section [3] will simplify our task significantly. 
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5.2 The proof of the lower bound in Theorem [5] 

We begin with step 1 of our heuristic, considering the size of Nj-if, £, Vj) for small r/. 
First we will need the following simple lemma. 

Lemma 24: 

For any t, (5 > and k > 0, 



Proof. Defining hs : 



l{^^ei-S,S)}d'^>k]<3e'/'-'/^^'\ 



Mby 
hsix) :-- 



\x\ if |x| > 5 

s J- ^ 
2 "•" 2<5 



^ + t if kl < (5 



we have, by approximating with C^ functions and applying Ito's formula, that 

6 



hsiit 



+ j h'sic,)d^s + Y^I 1 



{6e(-5,5)} 



ds 



(this function hs is often seen when studying local times of Brownian motion — for 
a full proof of the above see, for example. Lemma 4.11 of [27]). Also 



Thus 







- Jo h'si^s)d(s^ < ]p[e-/o'^^(6)d6-|/,>i(6)'d.]gt/2 < gt/2_ 



l|^,e{-5,5)}ds > A;) = P (hsiCt) -I- I h's{is)dis > ^ 



< 



\^t 



^ JO 
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k 



<P(I6I>^)+: 
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\^(6)d6>^ 



< 



J6I 



-fc/(45) _^ ]p 



- Jo h'siis)dis 



'k/{A5) 



< 3gt/2-fc/(4<5) 



establishing the result. 
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Lemma 25: 

For any continuous f with /(O) = and any e > 0, there exist rj>0, vyOjkyO 
and Ti such that for all T >Ti, 

n\NT{f,e/2,r,)\<uT)<e-^^. 

Proof. We first show that there exist r/ > 0, ki > and Ti such that 

]u G N{r]T) ■ u NTif,s/2,T])) < e'^^'^ VT > Ti. 
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Choose rj small enough that sup^grg^^i |/(s)| < e/4. Then, using the many-to-one 
lemma (at (*)) and standard properties of Brownian motion, 

P(3n e N{r]T) : u iVr(/,e/2,r/)) 

]u G Af(r/r), s < 7? : jX^sT) - T'if{s)\ > eT^l) 



< F{3u E N{riT) : sup |X„(s)| > sXyA) 

s<r]T 



< ^P(3n G A^(r/r) : sup |X„(s)| G [terV4, (A; + l)eTV4]) 






< V €""^^0 {{k+i)eT-^rdsp^ g^p i^^i g [;teTV4, (A; + l)erV4]) W 

k>i '^"f^ 



fe>i 

< j; -i= exp {{m/3ePrj - e^ {^2i^))kT^'i-^) 

for sufficiently small rj. For small r/ this is approximately 

Cexp ((m/3ePr/ - eV(32r/))T29-i) , 

which gives a decay of exp(— /jiT^''"^), which is more than the decay required. We 
now aim to show that for any ?? > 0, there exist v > Q and A;2 > such that 

¥{N{r]T) < uT) < e'^^'^ . 

Indeed, if we let n(t) be the number of births along the spine by time t, then certainly 

P(iV(r/T) < vT) 
< F{n{r]T) < uT) 






'0 

Lemma [24] shows that 



' ^{6e[-(4i//(/3r?))i/p,{4i//(/3r,))i/P]}ds > -riT \ < 3exp ( 



rjT r]T 



2 8(4i//(/3r/))Vp) 



so we have exponential decay in the first term provided that z/ < /3r//4P'^^; and since 
births along the spine occur at rate at least 'iv/r] outside the interval 

[-(4^/(/3r/))VP,(4i./(/3r/))i/P] 
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the second term is bounded above by the probabiHty that a Poisson random variable 
with mean 2z^T is less than vT. Let Y ~ Po(2i/T); then 

'^Y<vT = lexp(i^T)>cxp{y) < -^ 
SO 

P{Y < vT) < e'^'Ele-^] = ^uT+2.T{e.p{-i)-i) 

and this exponent is negative, so the second term also decays exponentially. Finally, 

P(|A^t(/, e/2, r/)| < uT) < P(3n G N{riT) : u Nrif, e/2, r?)) + ¥{N{r]T) < vT) 

and the proof is complete. D 

We now move on to step 2, using the results of Section |3] to bound the probability 
that we have a small number of particles strictly below 1. The bound given is 
extremely crude, and there is much room for manoeuvre in the proof, but any 
improvement would only add unnecessary detail. The proof of this lemma runs 
exactly as for the corresponding result in the p = Q case seen in [17j , with no extra 
technicalities; we include it again here for completeness. 

Lemma 26: 

/// G C^[0, 1] and K{f, s) > Vs G (0, 9], then for any e > and 6 > there exists 
To > and 7 < 1 such that 

P {\NT{f,e,d)\ < e^imT'^-'~ST^''~'\^ < 7 VT > Tq. 
Proof. Note that by Lemma [18] for small enough e > and large enough T, 

and hence 

P {\NT{f,e,e)\ < e^(/,e)T2'--5T^'-) < p (^zrier) < e-^^''-V2^ 

Suppose first that /'(O) = 0. Then K[Zt{9T)] = 1 and, again for small enough e. 



by Proposition 17 the set {Zt{9T),T > l,t G [1,0T]} is uniformly integrable. Thus 



we may choose K such that 



snpE[ZT{eT)l{z^^gT)>K}] < 1/4, 

T>1 



and then 



1 = E[Zt{9T)] = E[ZT{9T)l{z^^eT)<i/2}]+nZT{eT)l{y2<Zr{eT)<K}] 

+ E[ZT{eT)l{Zr{eT)>K}] 

< 1/2 + KF{Zt{9T) > 1/2) + 1/4 

so that 

¥{Zt{9T) > 1/2) > 1/(4K). 
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Hence for large enough T, 

p(|iVr(/,e,e)| <e^(/'^)^'^~'-^^'''') <1-1/{4K). 

This is true for all small e > 0; but increasing e only increases \Nj'{f,£,0)\ so the 
statement holds for all e > 0. Finally, if /'(O) / then choose g G C^[0, 0] such 
that 5(0) = g'{0) = 0, sup,<e |/ - ^l < e/2, K{g,(l)) > V,/. < ^ and K{g,e) > 
K{f, 9) — 5/2 (for small rj, the function 

(f\ _ / /(*) +at + ht^ + ct^ + dt^ if t £ [0, r]) 

^^*^'-\ fit) ifte[r?,i] 

will work for suitable a, b,c,d£M). Then 

F{\NTif,e,9)\ < e^'-f'^^)^''-'-'^"'-') < F{\NT{g,s/2,e)\ < e^MT"'-'-ST^^-'/2^ 

< 1 - 1/(4K) 

as required. D 

We are now ready to carry out step 3 of the heuristic. Again this runs exactly 
as in [T7] . 

Proposition 27: 

Suppose that f £ C^[0, 1] and K{f,s) > Vs G (0,0]. Then for lattice times Tj 
(recall that this means that there exists r > such that Tj — Tj^i = jr for all 
3 > V, 

limmf-^log\NT^{f,e,e)\ > K{f,e) 

almost surely. 

Proof. For a particle u, recall that u < v means that u is an ancestor of v and define 

NT{u,f,e,e) := {v G N{9T) : n < w, \X^,{tT) -T'^ f{t)\ < eT" Vt G [0,6]}, 
the set of descendants of u that are in Nj^if, £, 0). Then for 5 > and r] G [0, 9], 



T, 



«eAfT(/,£/2,r?) 

< n Fi\NT{g,el2,Q-ri)\ < e^(/.^)^^-^-^^^-^ 



riT 



u(^NT(f,e/2,rj) 

since, given J^r^T, {\^t{u, f,e,9)\ : u G A^t(/, e/2,r/)} are independent random vari- 
ables, and where g' : [0, 1] — t- M is any twice continuously differentiable extension of 
the function 

g: [0,9 -v] -^ K 

t ^ f^t + r])-f{r]). 
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If r] is small enough, then 



\Kif,e)-Kig,e-r^)\<6/2 



and 



K{g,s)>0 yse{0,9-r]]. 
Hence, applying Lemma [26} there exists 7 < 1 such that for all large T, 

< P {\NT{g,e/2,e- r?)| < e^(s,e-.)T^-^-5T2-V2J 



< 



7- 



Thus for large T, 



(41) 



Taking expectations in (41), and then applying Lemma 25, for small rj and some 



fjk > 0, for large T we have 



<P(|7VT(/,e/2,r/)|<i.r) + 7^^ 

The Borel-Cantelli lemma now tells us that for any lattice times Tj, j > 0, 

p[liminf-i3Tlog|iV,(/,£,e)|<K(/,e)-<5) =0, 

and taking a union over 6 > gives the result. 

We now move to continuous time using step 4 of our heuristic. 

Proposition 28: 

Suppose that f G C^[0, 1] and K{f, s) > Vs G (0,6*]. Then 

\unM-^\og\NT{f,e,e)\>K{f,9) 

almost surely. 

Proof. We claim first that for large enough j G N, provided that Ti < 1, 



D 



|iVT,,(/,e,^)|> inf |iVr(/,2e,^| 



(Zl^ve NT^{f,e,0),ue N{eT,+i) : v < u, sup |X„(0r) - X„(^T,)| > 



ei;^ 
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Indeed, ii v £ NTj{f,e,9), T G [Tj',Tj+i] and S G [0,6'T] then for any descendant u 
of V at time 9T, 

\Xu{S) - Tif{S/T)\ < \Xu{S) - Xu{S A eTj)\ + \Xu{S A eTj) - TU{^)\ 

+ \Tjfi^) - T]f{S/T)\ + \T]f{S/T) - T'if{S/T)\ 
<\Xu{S)-Xu{S^eTj)\+eT] 

+ T] sup |/(x)-/(y)|+ sup |/(x)||r/_^,-r/| 

x,ye[0,6l] xe[0,6l] 

k-?/l<l/r, 

< \X^{S) - Xu{S A ^r,)| + ^T] for large j; 

so that if any particle is in N^^ (/, £■, 0) but does not have a descendant in Nx{f, 2e, 6) 
then its descendants must satisfy 

sup \Xu{S)-X^{Tj)\>eTy2. 
This is enough to establish the claim, and we deduce via the many-to-one lemma 



plus Lemma 15 and standard properties of Brownian motion (see Proposition 4.15 



of [27] for a more detailed justification) that 



P(|iVT,(/,e,^)|> inf |A^t(/,2£,( 



< 4e™^^i ~ ™Pses(/,.) /o \9{s)\^As+ml3T^'i-\\f{9)\+(k+?.)e/2) -y ^-(keT 



,)V(80Ti) 



which, as in Lemma 25, is exponentially small in Tj. Thus the probabilities are 



summable and we may apply Borel-Cantelli to see that 



P(|iVT, (/, e,e)\> inf {NtW, 2e, 9)\ infinitely often) = 0. 



Now, 



liminf-2^1og|iVT(/,e,0)| <i^(/,0) 



< P I hminf ;^2^ log \Nt^ (/, 2e, 9)\ < K{f, i 



1-5-00 T 



. infTe[r„T,+i]I^^T(/,e,^)| \ 
hm mf , ^ , j — < 1 



which is zero by Proposition 27 and the above. D 



This gives us our desired lower bound for Theorem [5} 
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Corollary 29: 

For any open set A C C[0, 1] and 9 G [0, 1], we have 

\iui\ni-^\og\NT{A,e)\ > sup{K{f,t) : / G A,0o(/) > 0} 

almost surely. 

Proof. Clearly if sup{i^(/, t) : / G yl,0o(/) > ^} = -co then there is nothing to 
prove. Thus it suffices to consider the case when there exists f € A such that / G Hi 
and 6 < 9o{f). Since A is open, in this case we can in fact find f € A such that 
K{f, s) > Vs G (0, 6] (if K{f, (j)) = for some (p < 0, just choose rj small enough 
that (1 — r])f G A) and such that / is twice continuously differentiable on [0, 1] (the 
twice continuously differentiable functions are dense in C[0,1]). Thus necessarily 
sup{K{f, t) : f £ A, 6*0(7) > ^} > 0, and for any S > we may further assume that 
K(f,9) > sup{-fir(/, t) : / G A,6Q{f) > 6} — 5. Again since A is open, we may take 
£ such that B{f,e) C A; then clearly for any T 

NT{f,e,e)cNTiA,e) 



so by Proposition 27 we have 



hminf -^ log Nt{A, 6) > snp{K{f, t) : f e A, Ooif) > 9} - 6 
almost surely, and by taking a union over 5 > we may deduce the result. D 

5.3 The upper bound in Theorem [5] 

Our plan is as follows: we recall that we ruled out the possibility of any particles 



following unusual paths in Lemma 19 which allows us to restrict our attention to a 
compact set, and hence small balls about sensible paths. We then carry out the task 
of obtaining a bound along lattice times for balls about such paths in Proposition |30[ 



By expanding these balls slightly (using an argument similar to that in Proposition 



28 ) we may then bound the growth in continuous time; this is done in Lemma 31 and 



finally we draw this work together in Proposition 32 to give the bound in continuous 
time for any closed set D. 

For simplicity of notation, we break with convention by letting 



:= sup |/(s)| 
se[o,e] 

for / G C[0, 6*] or / G C[0, 1] (on this latter space, || • \\e is clearly not a norm, but this 
will not matter to us). We also extend the definition of N^iD, 9) to sets D C C[0, 9] 
in the obvious way, setting 

Nt{D, 9) ■= {u G N{9T) -.3/ gD with X„(tr) = T'^f{t) Vt G [0, 9]}. 

With a slight abuse of notation, for D C C[0, 1] and 9 G [0, 1] we define 

K{D,9):= sup K{f, 9). 
feD 
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We now attempt to establish an upper bound along lattice times for closed balls 



about functions outside F^- Recall the definition of F/v from Lemma 19 and that 



of i?Ar(e) from Proposition 23 



Proposition 30: 

Fix N £ N. For any closed ball D = B{f,e) C C[0, 1] about any f ^ Fjy, and any 
9 G [0, 1] and lattice times Tj, we have 

lim sup ^^ log \Nt^{D,9)\ < K{D, 9) + RN{e) 



rpz 



almost surely. 



Proof. Proposition 23 tells us that 



limsup^^logE[|iVT(Z),0)|] <K{D,9) + RN{e). 
Applying Markov's inequality, for any 5 > and p G [0, 2) we get 



r->oo 



^r 1 1 E[\Nt{D,9)\] 

< hmsup p^ log ^K^o,e)T^,-^^R^^e)T^,-.+ST^'^-^ ^ '^ 



so that for lattice times Ti , T2 , . . . we have 

00 

i=i 
and hence by the Borel-Cantelli lemma 

p(limsup^^log|Afr,p,0)| > K{D,9) + RN{e) + 6] =0. 

Taking a union over 5 > now gives the result. D 

We now check that an upper bound holds in continuous time. For e > and 
DcC[0,l], define 

D' := {/ G C[0, l]:3gGD with ||/ - g\\ < e}. 

Lemma 31: 

IfD = B{f,e) C C[0, 1] for some f Fn, then 

hmsup -^ log \Nt(,D, 9)\ < K{D', 9) + i?jv(2e) 
almost surely. 
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Proof. First note that for lattice times Ti,T2, . . ., 

\imsnp,^^log\NTiD,9)\ > K{D\e) + RN{2e) + 5\ 

<wUmsnv^^\og\NT^{D',e)\ >K{D',9) + RN{2e) 



j 



+ P lim sup -^ — T log sup —- — —^ — -T-r > 5 



] 



Clearly D^ = B{f,2e), so immediately by Proposition 30 



limsup^^log|7VT,(I?^0)| >K{D',e) + RN{2e) 
j 






and we may concentrate on the last term. We claim that for j large enough, provided 
that Ti < 1, for any T G [Tj,Tj^i] we have 

u G Nt{D, e)^3v<u with v G Nt^{D'', 9). 

Indeed, if n G iVT(D, 6) then for any 5 < OTj, 

\Xu{S) - T]f{S/T,)\ < \Xu{S) - ry {S/T) \ + \T]f {S/T,) - T'^f {S/T,) I 

+ TV{S/T,)-f{S/T)\ 

< r% + 11/11,(17+1 - TJ) + r^ sup |/(x) - f{y)\ 

x,ye[o,e] 
\x-y\<i/T, 

which is smaller than 2eT? for large j since / is absolutely continuous. 

We deduce that for large j every particle in Nt{D, 6) for any T G [Tj, Tj+i] has 
an ancestor in Nt {D'^,9). We now use this fact to ensure that Nt{D,6) cannot 
increase dramatically between times Tj and Tj+i. 

We temporarily need some more notation. For T > S > and u G N{S), let 
N{u, S, T) be the set of descendants of u born between times S and T. Also let 
Px be the translation of P under which we start with one particle at x rather than 
at the origin. Then, using the Markov property and the many-to-one lemma, for j 
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large enough, 



E 



sup {NriD, 
re[T„T,+i] 



J^eTi 



<E 






J^er, 



< Yl ^x^m) mien] 



u£Nt-{D 



.ml3 j;;^^ l^Pds 



< 



< 






\NT^{D^,e)\Ye 



m/39Ti(T/(||/|!e+2e)- 



A:>0 



A„-k^/(2eTi) 



since p < 2 this sum converges, giving 



E 



sup \Nt{D,( 



Tf 



9T, 



<\NTAD',e)\e^^^r^ 



where the 0(Tf'^) is deterministic. But pq = 2q — 2 and by Markov's inequahty 



sup r^ — TTT — -7- > exp oT- ^ 



<E 



E 



suPre[T,-,T,.+i] |^r(-D, 



J'f. 



er. 



\NTAD^,i 



exp(— (5r, 



2q-l-, 



Thus we may apply Borel-Cantelli to see that 



lim sup 2q-i log sup 



\NTiD,i 



>5\ =0. 



Again taking a union over 6 > gives the result. 



D 



We are now in a position to give an upper bound for any closed set D in contin- 
uous time. This upper bound is not quite what we asked for in Theorem [Sj but the 



final step — truncating K at 9q — will be carried out in Corollary 33 
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Proposition 32: 

If D C C[0, 1] is dosed, then for any 9 £ [0, 1] 

lhjisnp-^log\NT{D,0)\ < K{D,e) 

almost surely. 

Proof. Clearly (since our first particle starts from 0) Nt{D \ Co[0, 1], ^) = for all 
T, so we may assume without loss of generality that D C Co[0, 1]. Now fix 5 > 



and choose N (by Lemma 19) such that 



limsup log|JVT(-^Af,^)| = -oo, a.s. 



By Lemma 21 and the fact that Rn{2£) — t- as e — )• 0, we may choose e > such 



that K{D^,e) + RN{2e) < K{D,9) + 6. Then, by Lemma [20l for any N and some 
n (depending on N) and fk S C[0, 1] \ Fjv, k = 1,2, . . . ,n, 

\imsup-^log\NT{D,e)\ > KiD,9) + 6 
< P Aim sup ^^ log \Nt{Fn, e)\ > K{D, 6) + 5 



Y, IP (li^s^P ^4^ l°g l^r(^' ^' ^)l > ^(^'' ^) + «^(2e)) 



By our choice of N , the first term on the right-hand side is zero, and by Lemma 31 
all of the terms in the sum are also zero. As usual we take a union over (5 > to 
complete the proof. D 

Corollary 33: 

For any closed set D C C[0, 1] and 9 £ [0, 1], we have 

liuisnp-^log\NT{D,9)\ < sup{K{f,9) : / G D,9o{f) > 9} 

almost surely. 

Proof. Since \Nt{D,9)\ is integer valued, 

^log|iVT(A^)|<0 ^ ^log\NT{D,9)\ = -oo. 
Thus, by Proposition [30] if K{D,9) < then 

limsup ^^ log \Nt{D, e)\ > -oo j = 0. 

Further, clearly for <j) < 9 and any T > 0, if Nt{D,(P) = then necessarily 
Nt{D, 9) = 0. Thus if there exists (j) < 9 with K{D, 0) < 0, then 



limsup -2^ log |A^t(I),^)| >-oo) =0 



which completes the proof. D 
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Proof of Theorem\^ Combining Corollary |29| with Corollary |33| completes the proof. 

D 

Proof of Corollary^ Corollary ^ is proven in the same way as Corollary ^ (it is 
easily checked that the limsup and liminf of Theorem |4] agree on sets of the form 
D,,, := {/ e C[0, 1], 1/(1) -z\<e] and A,,, := {/ E C[0, 1], |/(1) - z\ < e] for 
e>0). D 

6 Proof of Theorem [8] 

There is much overlap in the proofs of Theorems |8] and [5} and so we will leave out 
many of the details and refer to earlier sections. Our task is also made easier by the 
fact that we are trying to prove a statement about one-dimensional distributions 
rather than the full sample paths of the BBM. 



Proof of upper hound in Theorem^ Lemma 19 tells us that for any A; > 0, we may 
choose A^ G N such that 

¥{Nt{Fn, 9)^(/})< exp {-kT^'^-^) (42) 

(recall that for large A'^, -F/v was a set of extreme paths that were very difficult for 



particles to follow). Now, in Proposition 23 we saw that for any / G Z? \ F/v and 
any(^e [0,1], 

limsup ^^logE[|iVT(/,e,<^)|] < K{D' ,^) + RN{e) 

where for each N, Rn{£) — )• as e — )• 0. 

Fix (5 > 0; choose k > such that —k < inf0<e K{D, (p) and N large enough 



that (42) is verified. By a similar argument to that in Lemma 21 (the upper- 
semicontinuity carries over easily to the function / i— )• inf0<0 K{f, 0)) we may then 
choose e small enough that 

inf K{D^, (j)) + RN{e) < inf K{D, (j)) + 6. 

(j><6 </'<^ 



Using compactness (see Lemma 20 ) , we can choose n G N and fi, f2, ■ ■ ■ , fn ^ D\ Pn 
such that 

DcFNUB{f,,e)U...UB{fn,e)cD'. 

Now for (j) < 6 and any set B, 

{Nt{B, </.) = 0} C {Nt{B, 9) = 0} 
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so, for T large enough (depending on the /i, . . . , /„) 

n 

F{Nt{D,9) ^ 0) < ¥{NT{FN,e) / 0) + ^P(iVT(/i,^) ^ 0) 

i=l 
n 

< P(iVr(Fw, ^) / 0) + 5Z i"f P(^t(/*, 0) / 0) 



i=l 



(/)<6» 



<e-'=^^^-^+5:infE[|iVH/, 



^-^ (b<e 



0<0 



Taking logarithms, dividing by T'^'^ "^ and letting (5^0 gives the desired result. D 

Proof of lower bound in Theorem [^ First note that it suffices to consider the case 
A = B{f,e) for / G C72[0,l]. Now 

f{Nt{a, e)^%) = F{ZT{eT) > o) 

1 



Zrier) 



> 



iT 



1 



}T[ZT{OT)\gT] 



As in the proof of Lemma 16 we can use Lemmas [14] and [15] to bound the spine 
decomposition; for any 6 > we may choose e small enough that 

Thus (since 5 > was arbitrary) it certainly suffices to show that 



lim inf Q^^ 

r-s>oo 






>0. 



But using the auxiliary random variables Ai,A2,--- and Si, S2, ■ ■ ■ from the proof 
of Proposition [17] we have 



T limsup Y. A^e-'fo- l^^l^ds < 00 > Q^ X]^"^"''"" < 



00 =1. 



T-s>oo 



We are now done by Fatou's lemma. 



\n=l 



D 
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7 Optimal paths: proofs of Theorems [4] and [7] 

We start with Theorem [7j We recah our optimisation problem. 

I / r«7./?ifc«^ip- ^ 

Kif,t):-- 




mf]\f{s)\P--f'{s[ 



ds if/eFi, 

otherwise, 
where the class of functions in which we would like to optimise is 

Hi := {/ : [0, 1] ^ M : f{t) = [ h{s) ds, h G L^[0, 1]}. 

Jo 

Optimisation Problem. For z € M, find 

sup K{f,l), 
feHuf{i)=z 

subject to 

K{f,e)>o ye € [0,1]. (43) 

By symmetry it is sufficient to consider z > 0. Furthermore, if 5 G Hi is any 



path satisfying the constraint (43), then \g\ E Hi also satisfies the constraint, has the 
same value at the end point, and K{g, 1) = -ftr(|5|, 1). So, without loss of generality, 
we can assume for existence that g > and drop the modulus sign from the definition 
of K. 

Optimal paths. For each s G [0, 1] define a path g G Hi as follows 

(i) For all t G [0, s), we set g{t) = r{t), the scaled position of the right-most 



particle (see (16)). Otherwise said, g is the solution to 

^g'itf = mf3g{t)P, g{0) = 0, t < s. 

(ii) For t £ [s,l], g satisfies 

g"{t)+m(3pgP-\t)=0, 
with initial condition specified by the fact that g and g' are continuous at s. 

We aim to show that, for each z G [0, z] the above optimisation problem has a 
unique solution given by the unique g^ constructed as above by picking the unique 
s = Sz £ [0, 1] such that 5^(1) = z. 



The strategy for our proof is that we will first show that gz defined above is the 

Apiece wise 



unique solution to the optimisation problem in the smaller class C'oiecewisc *- ^1 ^^ 



functions which are piecewise C^ on [0, 1] (note that the optimal functions gz have 
discontinuous second derivative at Sz)- We will do this by exhibiting a series of 
properties which a solution to the optimisation problem, if one exists, must satisfy. 
The unique function satisfying all of these properties will be our gz- We will then 

45 



show that Qz is also the optimum in the fuh class of functions Hi by showing that 
if there existed a better function in the larger class, there would also be a better 
function in our restricted class, and so obtain a contradiction. 

For notational convenience, in places where it will not cause confusion, we will 
write g for Qz- 

Lemma 34: 

Any optimal (in C^i^cewise) P"^^^ 9 is in fact C^. 

Proof. As g is C^if^cewise^ there exist = xq < xi < . . . < xk = 1 such that g is 
C^ on any interval {xi,XiJ^i),i = 0, . . . ,iC — 1. It is therefore enough to show that 
Vi = 1, . . . , -fC — 1, g'{xi—) = g'{xi+). To simplify notations we just write x for Xj in 
the following. 

Suppose that g'{x+) < g'{x—). We will show that it is possible to construct a 
better function. Choose e > 0, and consider the function g defined by taking g, and 
interpolating linearly on the interval [x — £,x + e]: 

'git), iG[0,x-e) 

g{t) := < g{x - e) + ^^^^^ {g{x + s) - g{x - s)j , tG[x-e,x + e) 
^g{t), tG[x + £,!]. 

For t £ (x — e,x), we have g'{t) = g'{x—) + 0{e); similarly, for t € (x,x + e), 
we have g'{t) = g'{x+) + 0(e). To simplify notation, we define g'_ := g'{x—) and 
g'j^ := g'{x+). Using Taylor expansions again, we have 

K{g, 1) - K{g, 1) = £^' l^-g'{tf - \g' [tf^ dt + 0{e^) 
= \{g'_-g'^f + 0{e')>Q 



for all sufficiently small e > 0. A very similar calculation shows that g satisfies (43) 
for all t G (x — e, X + e), which proves the result. The case g'j^ < g'_ is settled in the 
same manner. D 

Lemma 35: 

// g is an optimal trajectory then, on any interval I of [0, 1] such that K{g, t) > 
for all t € I, one has 

mppgP-^ + g" = ^. (44) 

Proof. Starting from an optimal trajectory g, consider the deformed trajectory g+eh 
where e is small and /i is a sufficiently smooth function with /i(0) = h{\) = so 
that in particular the end-point is still fixed at z. Using a Taylor expansion and an 
integration by parts we have 



K{g + eh,t) = K{g,t)+ I [ml3{g + shy - mpg^ - g'h' e\ds + 0{e^), 

Jo 

= K{g,t)+e I [mppgP-^ + g"]hds - eg' {t)h{t) + 0{£^). 
Jo 



(45) 
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Assume that K(g, t) > on some interval /. Then for any ti < t2 in /, there exists 



c > such that K{g,t) > c for all t € [^1,^2]- Assuming that (44) does not hold, 
choose h to be of the same sign of m/3pgP~^ + g" in t G [^1,^2] and zero everywhere 
else. Then, for e small enough, K{g + eh, 1) > K{g, 1) and K{g + eh, t) > for all 
t (one simply needs to choose e so that eg'{t)h{t) — O(e^) < c for all t E [ti,t2]-) 



Therefore, if g does not satisfy (44), then g + e/i is a better path. D 



Recall from (16) that 



1 



r{s) = -^—{2-p) 



describes the limiting shape of the boundary of the trace of the rescaled BBM. It 
solves 

-r'{sf = ml3r{sY, r(0) = 0. 

We are now going to prove that it is not possible to find a path / along which the 



condition (43) remains valid and f{s) > r{s) for some s G (0, 1). Indeed this would 



clearly yield a contradiction with Theorem [T} 
Lemma 36: 

-^2 

'' piecewisC 



Let f e Cp,ece™.e' /(O) = 0, / > 6e such that ^ holds. Then 



f{s)<r{s),ysG [0,1]. 



Proof. Suppose that there exists / S C'piecewise' /(O) = 0, / > such that (43) holds 



and such that /(sq) > ^(so) for some sq G (0, 1). We construct g € C'piecewise' 5(0) — 
0,g > as follows: let si = so if K{f, so) > and otherwise pick si > so so that 
9{so) > r(si). 

- On [0, So] take g = f, 

- On [so,si] take g{s) = /(so), 

- On [si, 1] g is the solution of 

lg'{sf = mf3gisy, g{si) = f{so). 

Observe that g has the properties that g{l) > r{l), K{g, si) > and K{g,l) = 
K{g, si). According to Theorem Isl for any e > 

2+P 

liminf T 2-p log \NT{B{g, e),l)\> K{g, 1) = K{g, si) > 

T-)-oo 

which contradicts Theorem [TJ D 
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Now to prove the first line of Theorem [7j it is sufficient to show that 

sup{K(/,l),0 < f{s) < r{s),f{l) = z}< K,,,,{z). (46) 

Take / such that f < r and 3t < 1, K{f, t) < 0. Then / cannot be optimal among 
paths which stay below r: one can easily construct a better path / staying below r 
by choosing / = r up to r~^(/(t)), constant between that point and t, and equal to 
/ thereafter. 

Lemma 37: 

Define Sz{g) := inf{s G [0, 1] : g{s) / ^'(s)}. Then any optimal path g solves 

ml3pgP-^ +g" = 
on the interval s G {sz{g), !)■ Furthermore, K{g, s) > for all s G {sz{g), 1). 

Proof. If Sz{g) = 1, there is nothing to prove, so we assume Sz{g) G [0, 1). The first 
thing to observe is that any point (t,y) G {f G [0, 1],0 < y < r{t)} can be reached 
by a path / such that f(t) = y and K{f, t) > 0. If y > just let / = r until 
r(s) = y and then let / stay constant (K stays at until s and then accumulates 
some positive growth until t). If y = it is equally easy to check that there is a 



path which reaches this point, fulfills (43) and has a strictly positive K at the end. 



Suppose now that there exists t > Sz such that g{t) = r{t). Then there must 
exist [a, b] C [s^, t], a < b, such that g{a) — r{a) = g{b) — r{b) = and r{t) — g{t) > 
for all t G (a, 6). But by the previous observation this means that K[g, t) > for all 



t G (a, b) and therefore by Lemma 35 g must be strictly concave on (a, b). But since 



R is strictly convex this is a contradiction. Thus we must have K > Q on {sz, 1] and 



we conclude by Lemma 35 again. D 



Lemmas 35 , 34] and [37] show that any positive solution gz to the optimisation 
problem (in C^iecewise) ^^ such that there exists Sz such that 

(i) gz=r on [OjSz], 

(ii) gz solves g" + mfipgl" = and gz{l) = z 

(iii) g'z is continuous at Sz- 

Therefore we conclude that the unique positive solution to the optimisation problem 
in Cpjg(,g^jgg is gz- (Note that it is easily seen that Sz > for p > 0.) 



Lemma 38: 

irfini y 'if n -/« nn nnHmnl rinfh rn)PT C , . 

piecewise 



For any endpoint z, if g is an optimal path over C^jecewise ^^^'^ 9 ^-^ ^^^^ optimal over 



Hi. 

Proof. We have shown the existence of an optimal path gz over all paths in Cpig^g^jg^; 
suppose there is a better path f £ Hi. That is, /(O) = 0, /(I) = z, f > 0, 
K{f,t) > Vi G [0,1] and K(/,l) > K{gz,l) + e for some e > 0. For each 
n G N, define /i„ : [0, 1] — ^ M by setting hn{k/n) = f{k/n) for ah A; = 0, 1, 2, . . . , n 
and interpolating linearly elsewhere. Then each /i„ is piecewise linear and hence 
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certainly piecewise C^; and since hn agrees with / at each k/n and hnear functions 
minimise derivatives, we have 



2 

n n' 



/ K{sfds< [ fisfds vtejo,-, 

Jo Jo in 

Now by choosing n large we may insist that 

[ h'^isfds < I f'isfds + e Vt G [0, 1] 
Jo Jo 

and 

/ K{srds> [ fisyds--^ VtG[0,l]; 
Jo Jo 2m/3 

but then /i„ is a function in Cpjg^g^jgg satisfying all the required properties and with 
K(hn,l) > K{f,l) — e > K{gz,l). This contradicts the assumption that gz was 
optimal in C^.^.^^i^g. D 

All that is left to prove the first part of Theorem [7] is to show uniqueness in Hi . 

Proof of uniqueness. Now fix z > and suppose that 3h E Hi such that /i(l) = z 
and K(h, 1) = K{gz, 1) and OQ{h) > 1 but that h ^ gz- Take s S [0, 1] such that 
h{s) 7^ gz{s)- By rescaling time by 1/s and considering the endpoint h{s), by Lemma 



38 we can find some positive piecewise C^ function /i ending at /i(s) with growth 
rate K{f,s) > K{h,s). Equally, there is an optimal positive piecewise C^ function 
/2 amongst functions beginning at h(s) and ending at z; this is not immediate 
from our results as we have not considered starting from anywhere other than the 
origin, but our proofs easily carry over with no extra work. Since these two optimal 
growth rates (from to h{s), and from h{s) to z) are achieved by positive piecewise 
C^ functions, there is a positive piecewise C^ function / such that f{s) = h{s) 
and K{f, 1) > K{h, 1) = K{gz, 1). This contradicts the uniqueness of gz amongst 
positive piecewise C^ functions. 

Note that for z = and p > 0, the solution is not unique: the positive function 
go and the negative function —go are both optimal. D 

We now turn to the second part of Theorem [7] which concerns the total population 
size. 

Lemma 39: 

There exists a unique Za.s. > such that 

K^.s. ■■= i^a.s.(^a.s.) = sup K,.sXz) = SUp {K{f, 1),/ € C[0, liOoif) = Oo}. 

z 

The total population size satisfies 

lim — 2+f^og\N{T)\=Kii,s, almost surely, 

T-5-00 2^2^ 
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where one finds 



0, 



3p-2 



+ 



V2rnP 



and 



K, 



2-p 



2 
2-p 



2+p 



m(3zl^^ 



Proof. By Lemma 38 we may assume without loss of generality that g £ C'piecewise- 
Let e > be small and g be any such function satisfying (43). For f G [1 — e, 1] we 
have g{t) = g{l) + 0{e), and g'{t) = g'{l) + 0{e). Therefore 

^' (rnPgity - ^^'(t)') dt = e{m(3giir - ^^'(l)') + 0{e^). (47) 

'g{t), tG[0,l-e) 



Now consider the function 



g{l-e), te[l-e,l]. 



Note that g satisfies (43) because g does, and 



K{~g, 1) = K{g, 1 - e) + em/3g(l - ef 

= K{g,l-e) + empg{ir + 0{e^) 



(48) 



But comparing (47) and (48), we see that K{g, 1) > K{g, 1) for all sufficiently small 
£, unless g'{l) = 0, hence the first part of the result. 



For s > s, one has 



^giisf + mf3g,{sr = mf3A 



where A is some constant. For z = Za.s. one has g'~ ^ 0, g'~ (1) = and A 



Therefore (for s > Sz), 



9k M 



y^2ml3 



(49) 



(50) 



V^ls. - Qz,,,, (s)P 

We integrate the above expression with respect to s from Sz^^ to 1 and we make 
the change of variable x = gz^^ (s)/^a.s.- We obtain 



?l-p/2 



a.s. * 

1 



dx 



2-i/p \/l - XP 



V^i^ 



Sz, 



(51) 



For the lower limit on the integral, we used that 29'zi^z) = innl3gz{szY ^^"^ hence 
9z{szY = A/2. For z = Za..s. this gives gzisz^,^,) = Za.s.2~^/P. The value of Sz^^ then 
comes from the explicit expression gzis) = r(s) when s < Sz^^ . 

The expression for K^^s. comes from Theorem |9J which will be proved in Section 
8. D 
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The solution to the unconstrained optimisation problem given in Theorem |4] is 
now simple in light of the work above. 



Proof of Theorem^ Inspecting the proof of Lemma 35, we see that the solution to 
this unconstrained problem amongst functions which are piecewise C^ is given by 
the function hz satisfying 

h"{u) + m(3pW-^ = 0, 
/i(0) = 0, 
h{l) = z. 



A similar argument to that in Lemma 38 then shows that this function is also optimal 



over Hi, and another similar to that in Lemma 39 gives that the optimal z is that 
with h'^{l) = 0. 



The value of % is obtained in the same way as Za.s. except that (50) must be 
integrated from to 1 rather than from Sz to 1. The value of K^ also comes from 
Theorem [S D 



As an aside, we note that 



i^E(0) = 2 2-pi^^. (52) 



This can be seen by remarking that the optimal path /iq is symmetrical around 
s = 1/2 and hence that h'Q{l/2) = 0. The trajectory h^ up to s = 1/2 is therefore 
the trajectory maximising the total population at time T/2 and, given the total 

_2+p ^ 

population growth rate, K^{^) = K{hQ, 1) = 2 x K{hQ, 1/2) = 2x2 '^-p K^ which 



is the same as (52). 



8 Further properties of the optimal paths 

In this section we will prove Theorem ^ which states that both growth rate Ke{ 



and Ka.s.i^) ^■re solution of the same differential equation (22). We will give only 
one demonstration for both quantities, highlighting where necessary the differences 
between the two cases. We write in a generic way K{z) for either quantity K^{z) 
and i^a.s.(-z)- Similarly, fzis) stands in this section for either optimal path gz{s) or 
hz{s) defined respectively in Theorems |4] and [7J One can write 



K{z) = -^ / f'z{sfds + mP I fzisfds. (53) 

J Sz " Sz 

where fz{s) is a solution of 

f':{s) + mMzisf-^ = 0, /.(I) = z (54) 



and 



Sz = 0, fz{sz) = in the expectation case, 

(55) 
fz{sz) = t{sz)-, f'zi^z) = f {sz) in the almost sure case. 
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(In the almost sure case, the optimal path is equal to r{s) for < s < s^ where Sz is 
a unknown quantity which has to be solved for. Recall that r{s) is the trajectory of 
the almost sure rightmost particle; on this trajectory the population does not grow, 



which is why one can start the integrals in (53) from Sz- 



We begin with a simple lemma showing monotonicity of the optimal paths in z. 

Lemma 40: 

The optimal paths and their derivatives are monotone in z. That is, if < w < z, 

Us)<Ms) and fUs)<f'zis). 

Proof. Suppose there exists uq such that /«,(no) > fz{uQ). By the intermediate 
value theorem there exists s £ {szV Sw,uo) and t G {uq, 1) such that /«,(s) = fz{s) 
and fwit) = fz{t)- Then by our characterisation of the optimal functions in terms of 
solutions to differential equations, we must have /«,(u) = /^(u) for all u € [s, t\. This 
contradicts the fact that /«,(tio) > fz{uo). A similar proof works for /^ by considering 
hypothetical points ui < U2 such that fw{u2) — fwi^i) > fziu2) — fw{ui)- D 



Lemma 41: 

One has K'{z) 



-/^(l)- 



Proof. By Lemma 40 we may differentiate (53) with respect to z. One gets 
K'{z) 



f'zis)^is)ds + m^p 



fz{sr'^{s)ds 



dsz 
' dz 



1 



mPfziszT - -f'zisz) 



(56) 



The third term in the right-hand side is zero, either because dsz/dz = (expectation 



case) or because the square bracket is zero (almost sure case, see ( [19| ) with (55)). 
Integrating the first term by parts leads to 



K'{z) 



[f':{s) + m(3pf.{sr-'] ^{s)ds-f:{l)^{l) + fi{sz)^{sz). (57) 



dz 



dz 



dz 



The integral is null because of (54). As fzi^) = z for all z, one has dfz/dz{l) = 1 in 



the second term. The third term is also null because dfz/dz{sz) = 0. This is trivial 
in the expectation case as fz{sz) = /^(O) = and it is also true in the almost sure 
case becausqj/2(s) is independent of z up to s = Sz. D 



tegrating the second integral by parts and applying Lemma 41 , one gets 



Proof of Theorem^ We now establish relations between the integrals in (53). In 
.e second inte 

f f'M^'^s 

J s. 



-zK'{z)-fz{Sz)f'z{Sz) 



fz{s)f':{s)ds. 



(58) 



^To be more verbose, as fz(sz) = i'{sz) one has 



dz 



■r'{sz) 



Az 



oz dz 



Using f'{sz) — r'{sz) gives d}z/dz{sz) — 0. 
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Then, using (54) 



/ fUsfds-mf3p f Msrds = -zK'{z)-Us,)fUs.). 

Jsz Jsz 



(59) 



We now multiply (54) by /^(s) and integrate: 



^f,{sf + m/3Msr = c. (60) 

The integration constant c can be obtained by evaluating at s = 1 or at s = Sz: 

c = '^K'izf + m/3zP = \as.? + m^Us^y. (61) 



Then, integrating (60) between Sz and 1, 



^ / fiisfds + mpf Msrds = il-sz)c 



'-K\zf + mpz'' 



-f'^{szf + m(3Uszr 



■ (62) 



where both expressions for c were used. Now, the right linear combination of (59) 



adding, and using (53), one gets 



and (62) gives K{z): multiplying (59) by -2/(2 +p) and (62) by {2-p)/{2 + p), 



K{z) 



2zK'{z) + {2-p) [K'{zf/2 + mPzP] 



+ 



2 + p 

2Usz)asz) 



{2-p)sz[aszf/2 + mMz{szY] 



2 + p 



■ (63) 



The second term in the right-hand side is zero. This is trivial in the expectation case 
as Sx = and f{sz) = 0, and can easily be shown in the almost sure case: replace 
all fz by r using ( |55[ ), replace ml3r{szY by another r'{sz)'^ /2 from ( |19[ ) and after 
simplification one gets that the term is zero if r'{sz)/r{sz) = 2/[(2 — p)sz], which 



is true as can be seen from differentiating the logarithm of (16) with respect to s. 



Then, one checks easily that (63) is equivalent to (22). 



Note that for the optimal endpoint z one has K'[z) = and hence obtains 

2-p 



K{z) 



2+p 



mfSz^, 



(64) 



as stated in Theorems |4] and [7l 



It now remains to prove that (22) can be rewritten, for z > 0, as in (^27^ with a 



plus sign in front of the square root. For this, it is sufficient to show that K'{0'^) > 0. 
For p > 0, we have seen both in the expectation and the almost-sure case that 
the non- negative optimal path /o(s) going to the origin is not identically zero but 
goes away from the origin to take advantage of higher branching rates. Consider this 
positive optimal path and let Zm = maxg /o(s) > 0. For < z < Zm, consider the 
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path / which is equal to /o up to the point where z is reached for the last time and 
which is identically equal to z past that point. Clearly, one has i^(/, 1) > K{fQ, 1) 
and 6*0 (/) > Ooifo)- This implies that K{z) > K{0) for all positive z smaller than Zm 
and hence that K'{0~^) > 0. (Note that this argument holds both in the expectation 
and almost-sure cases, although the value of Zm is not the same.) For p = 0, optimal 
paths in expectation and almost-sure cases coincide, one has the explicit solution 
K{z) = mfi — z"^ /2 and K'iQ) = 0, the square root in (27) is zero and the sign is of 



no importance. D 
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